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Abstract 


Part I of this report describes a method for tidal analysis based on 
high and low water observations and also discusses the results of various 
test runs. Part II gives detailed information on certain aspects of the 
calculations and Part III is a user's manual for the relevant computer 
program. 


Wherever possible, terminology, computations and computer input and 
output formats are in accord with two earlier reports (Foreman, 1977, 1978) 
dealing with equally spaced tidal height and current data. 


Users who wish to receive updates of the computer programs should send 
their names, addresses and type of computer used for implementation to the 
authors. 
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1. INTRODUCTION. 


It is customary nowadays, in tidal heights analysis, to use hourly values 
of surface elevation obtained either from a digital tide gauge or by sampling 
a continuous record of water level. However, there is a mass of historical 
tidal data in the form of timesand magnitudes of successive high and low 
water levels, and even today there are situations in which it is convenient 
to record tidal elevations in this form. The fact that the time intervals 
between successive tidal extrema vary considerably means that most standard 
computer programs (e.g. Foreman, 1977) used for tidal (harmonic) analysis of 
conventional hourly-sampled data are not applicable. Yet these are generally 
based on least squares fitting of tidal frequency constituents to the data 
and there is no reason, in principle, why irregularly-spaced data cannot be 
analysed in the same way. In such an approach to analysis of high and low 
water levels which, for convenience, we term "high-low analysis", there is 
an upper limit to the tidal constituent frequencies which can be included. 
Basic sampling theory shows that constituents which are sampled fewer than 
two times per cycle become irretrievably confused with lower frequency 
constituents - a type of error known as "aliasing". Knowing the times and 
magnitudes of low and high waters should be somewhat better than simply 
having two elevation samples per tidal cycle, since the time derivative of 
surface elevation is zero at extrema. Assuming that, on average, four extrema 
occur per (lunar) day, it may be possible, using all the information implicit 
in the extremal amplitudes and times, to resolve constituents with periods as 
low as six hours. A conservative position is taken in the test cases 
described later in §3, where no attempt is made to determine constituents 
above semi-diurnal frequency. It is shown, however, that failure to use the 
derivative information leads to less accurate estimation of semi-diurnal 
constituents - a result probably attributable to aliasing of terdiurnal and 
quarter-diurnal constituents. This effect will be more noticeable in 
locations where higher-frequency constituents are relatively large, for 
instance, in shallow water. 


Although high and low water observations have been taken over many months 
at a few locations, a record length of one month was taken in the numerical 
tests as being more typical. Consequently, low-frequency constituents (e.g. 
fortnightly, monthly) were omitted from the analysis, since they often show 
considerable variation from month to month due to meteorological effects; 
there is no theoretical obstacle to their inclusion, however. 


Visual observations of tidal extrema often cover daylight hours only, in 
which case either one or two observations per day are missing and, at some 
sites, no observations at all were made on Sundays. These defects affect the 
accuracy of high-low analysis, and use of longer records is the obvious 
countermeasure. The missing data points cause no practical difficulty in 
least squares fitting procedures originally set up for unevenly spaced data. 


A most important consideration in analyses based on high and low water 
data is the reliability of the observed times. Ideally, observations of water 
level should be made from some time before the extreme value until some time 


afterward, and then plotted so that the peak value and its time of occurrence 
can be estimated accurately. But, even so, short period waves introduce error 
into all the observations and it is reasonable to expect a minimum random 
error in times of high and low water of the order of a few minutes; the 
effect of such random timing errors is discussed later. If the correct 
observation procedure, outlined above, is not followed, the times of extrema 
made by most observers tend to be several minutes late. Consistent errors of 
this type affect the phases of constituents by a calculable amount, since the 
period of each constituent is known, and corrections are possible if some 
estimate of the timing error can be made. 


Previous work on tidal analysis of high and low water observations ranges 
from the traditional sequential methods discussed by Schureman, 1958, which 
yield only moderately accurate estimates of the principal constituents, to 
least squares analysis (e.g. Zetler, Schuldt, Whipple and Hicks, 1965) 
essentially similar to that used here. The purpose of the present report is 
to combine an efficient least squares algorithm with the full suite of nodal 
modulation, astronomical argument correction and full inference calculations, 
while maintaining maximum compatibility of terminology and format with other 
tidal analysis programs currently in use at the Institute of Ocean Sciences. 


2. METHOD OF ANALYSIS. 
2.1 Choice of Constituents. 


The magnitude of tidal constituents usually does not vary rapidly with 
latitude or longitude and, if no. prior information is available from the site 
in question, it is fairly safe to assume that the same constituents will 
dominate the tidal behaviour at the new location as at other known stations in 
the same geographical area. In the numerical tests below, selection of first 
magnitude constituents was simple, since conventional harmonic analysis 
results for the test site were already available. The diurnals 0, and K, and 
semi-diurnals No, My and S. were, clearly, most significant; in certain 
tests, the diurnal P, and semi-diurnal K> were also included, being inferred 
from K,; and S»2, respectively (see 85.3). Constituents outside the diurnal to 
semi-diurnal frequency range were excluded for reasons discussed in 81. A 
constant term, Z,, was included as a matter of course, since surface 
elevations are seldom measured directly relative to mean water level. It may 
be noted that, with a record length of one month, the system of equations to 
be solved is substantially over-determined and there is no difficulty in 
adding minor constituents if desired. 


2.2 Least Squares Fitting of Sinusoidal Constituents to High and 
Low Water Data. 


Assuming that a sequence of high and low water observations, Y;> and the 
corresponding times, t., 1 = 1,...N, at which they occurred are given, we wish 
to find a function 


M 
t ech .* Z te ach . 
y(t) . 2, A, cos an(o,t 3) (1) 


in which the constituent frequencies o. and the number of constituents M are 
specified beforehand, but the amplitudes A. and phases F remain to be chosen 
so that the values y(t.) of the fitting function at the sampling instants t, 
agree as well as possible with the contemporaneous observed elevations Yas 
that is, 


M 
Y; -[A, + 4, a cos alot, me +) Fhe eas Oh | ip Toexen (2) 


Further, at the observation times t,, the time derivative y'(t) of the 
fitting function should be approximately zero, that is 
M 

y (t.) = - ve 2n0 5 A. Sin an(o st, - $5) = F « 0 (3) 
The fitting errors Eas 53s cannot be reduced exactly to zero when the 
number of arbitrary constants (2M + 1) in the expression for y(t) is less 
than 2N, the number of equations (2) and (3) to be satisfied. A commonly 
adopted compromise in such over-determined problems is to minimize the sum 


of the squares of errors at the observation times, which means, in the 
present case, choosing the A. and ¢. so as to minimize the error function 


E = 2 ifvep] [wr cep] x {., + vse} (4) 


that is, to find a "least squares fit" to the available data. The inclusion 
of an arbitrary positive weighting coefficient w in (4) permits control of 
the emphasis to be placed on satisfying the zero derivative condition 
compared to that placed on having y(t) fit the observed elevations accurately. 
For instance, w = 1.0 indicates that equal emphasis is given to both 
conditions, whereas w = QO means that the requirement that y'(t) should be 
approximately zero at each t. is simply ignored. 


Details of the algorithm used for numerical minimization of E in (4) are 
given in 84. 


3. NUMERICAL TESTS. 


To test the effectiveness and accuracy of high-low analysis by least 
squares fit, some numerical experiments were carried out on surface level 


observations collected in 1974 at Prince Rupert, B. C., a fairly typical 

West Coast port with relatively deep approaches, where the tide is pre- 
dominantly semi-diurnal with significant diurnal contributions. In order to 
have some basis for judging the high-low analysis, four conventional harmonic 
analyses of hourly heights were carried out first. 


3.1 Harmonic Analyses. 


Analysis 1. 


A full 12-month hourly height harmonic analysis (Foreman, 1977) for 1974 
is listed in Appendix 1. The constant component plus the major diurnal 
and semi-diurnal constituents from this 68-constituent analysis form the 
first row of Table 1. That there is little non-tidal contribution to 
water level variation at Prince Rupert is evident from the fact that the 
residual elevation, after removal of tidal constituents found in Analysis 
1, had an r.m.s. value of 0.13 m, which is approximately 2% of the tidal 
range. 


Since seasonal variation at Prince Rupert is very slight and one-month 
analyses were the principal topic of interest, all subsequent tests were 
confined to a single arbitrarily-selected month - January, 1974. 


Analysis 2. 


Hourly heights for January, 1974 were analysed for the suite of 36 
constituents which can be resolved from a 31-day record using a Rayleigh 
(resolution) criterion value of 0.97 (Foreman, 1977, p.16). Constituents 
Pi, ©3s Vo and Ky, which can only be obtained by inference in a one-month 
analysis (see §5.3), were omitted. 


Analysis 3. 


Though otherwise similar to Analysis 2 above, this case used 40 
constituents, P;, 91, vo and Ky being inferred from K,, Qi, No and So, 
uSing inference constants calculated from the one-year Analysis 1. 


Comparison of Analyses 1, 2 and 3 shows clearly how much more accurately 
the constituents K,, No and S» can be estimated from monthly records wiien 
' inference is used. Of course, in this example, the inference constants are 
optimum, since they were calculated from a longer record at the same site; 
in practice, inference constants may have to be estimated from data at a 
neighbouring site, which can result in a less marked improvement. 


Since tidal constituent frequencies are not harmonics of a single 
fundamental frequency, the results obtained by least squares fit analysis 
depend on the number of constituents included. The effect on the major 
constituents of omitting many of the minor constituents can be seen on 
comparing the 40-constituent Analysis 3 with the following : 
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Analysis 4. 


This test is a monthly harmonic analysis identical to Analysis 3 but 
for the fact that only five major, plus three inferred constituents and 
a constant term are included in the least square fit. These are the 
constituents included in most of the high-low analyses, described below. 


3.2 High-Low Analyses. 


The heights and times of high and low waters used in the following tests 
were taken from a strip-chart record. The times are estimated to lie within 
3 minutes of actual high or low water; heights are correct to within +1.5 cm 
(0.05 ft). As explained in §2.2, the relative importance accorded to the fact 
that the time derivative of elevation is zero at high and low water is 


reflected in the magnitude of the weighting coefficient in the leas t-Square- 
fit procedure. 


Analysis 5. 


This is a least square fit of five major constituents to 119 consecutive 
high and low waters (times and magnitudes) at Prince Rupert in January, 


1974. No inference of Pj, vy and Ky was carried out. The zero derivative 


information was accorded the same significance as the elevation readings, 
i.e. the weighting coefficient, w, was taken as unity. 


Comparing Analysis 5 with 2, it is clear that high-low analysis is 
capable of determining constant, diurnal and semi-diurnal constituents with 
very satisfactory accuracy. That inference is a useful additional feature 
in high-low analysis is verified in the following test : 


Analysis 6. 


This high-low analysis of the data already examined in Analysis 5 
employs the same five major constituents, but also infers Pisvo and Ko, 
using the same inference constants as in Analyses 3 and 4. 


Comparing these results with Analysis 5, and referring also to the 
constituents found in the one-year Analysis 1, one finds the same degree of 
improvement in constituents K,;, No and S» as when inference was introduced 
into the harmonic analysis (Analysis 3 versus Analysis 2). It can be 
concluded that inference is worthwhile in high-low analysis whenever 
reasonably reliable inference constants are available. 


Analysis 7. 


This test is identical to Analysis 6, except that the zero derivative 
information is not used (zero weighting coefficient) 
as in the foregoing test. 


. Inference was applied 


The results of Analysis 7 obviously differ from Analysis 3 rather more 
than Analysis 6 does; in other words, omission of derivative information 
impairs high-low analysis somewhat. Nevertheless, Analysis 7 is close enough 
to Analysis 3 to be considered quite satisfactory. This is an important 
conclusion, since Analysis 7 is probably representative of the results which 
can be expected from irregularly-spaced sets of elevation observations which 
are not necessarily extrema. In fact, the computer program used for high-low 
analyses can be used unchanged for any arbitrarily-chosen set of observations, 
provided the derivative weighting coefficient, w, iS zero. 


3.2.1 High-Low Sequences with Gaps. 


High and low waters which occurred during hours of darkness are often 
missing from records taken visually. In some cases, no readings were made on 
Sundays. In order to find out the effects of such gaps on the accuracy of 
high-low analyses, the data already used in Analyses 5 through 7 were 
modified by deleting observations made between 9 p.m. and 5 a.m. (local time) 
on weekdays and Saturdays, and all day Sunday. This reduced the total number 
of observations for January, 1974 from 119 to 73. 


Analysis 8. 


The reduced set of observations described above was subjected to high- 
low analysis for 5 major and 3 inferred constituents, with unit weighting 
for the derivative information. 


The results in Table 1 agree slightly less well with Analysis 3 than the 
various earlier analyses based on 119 observations but, nevertheless, are 
remarkably close. 


3.2.2 High-Low Analysis with Timing Errors. 


In order to simulate the effects of errors in timing when high and low 
water are estimated purely by eye, the genuine data already used in Analyses 
5 to 8 were altered by adding randomly-generated timing errors, uniformly 
distributed in the range -15 to +15 minutes, to the individual observation 
times. The magnitudes of the observed high and low waters were not altered. 


Analysis 9. 


A high-low analysis of the data with simulated timing errors, as 
described above, was carried out using inference and with a weighting 
coefficient of 1.0 for the zero derivative condition. 


Analysis 10. 
The preceding analysis was repeated with a weighting coefficient of 


zero, i.e. the fact that the observed elevations were known to be 
maxima or minima was ignored. 


The results of Analyses 9 and 10 show that the former gives estimates 
closer to the true tidal content of the data as defined by Analyses 3 and 4. 


3.3 CONCLUSIONS. 


It is clear from the numerical tests in §3 that "high-low analysis’, 
consisting of a least-squares fit of major diurnal and semi-diurnal tidal 
constituents to high and low water levels can yield very satisfactory 
estimates of amplitudes and phases of the constituents involved, at least 
for records about one month in length. It is expected that minor constituents 
of semi-diurnal and lower frequencies can be resolved where longer records 
are available. However, higher frequency constituents are best considered as 
noise, since even taking the derivative information into account, the upper 
limit placed by sampling theory on the resolvable frequencies iS somewhere 
between semi-diurnal and quarter-diurnal. 


The tests analyses indicate that 


(i) it is always best to use rather than to ignore the information 
that the observations of elevation are also extrema, 


(j7) the estimates of some major constituents are improved by 
"inferring" certain subsidiary constituents, 


(iii) a lack of night-time observations and occasional missing days 
impair high-low analysis only slightly, and 


(iv) high-low analysis is fairly insensitive to errors of several 
minutes in the observed times of high or low water. 


PART II. DETAILS OF PROCEDURES. 


4. Least Squares Fit with Modified Gram-Schmidt Algorithm. 


The original system of overdetermined equations which gives rise to the 
least squares fit problem is 


M 
y(t.) Tent A, + 2 A, COS an(o.t, - ;) cheer 0 
M 1 = TAR N (5) 
y (t.) Se & eno 5A. sinen(o.t, - o;) = 0 


in the notation of §2. It is convenient to change variables to Co = Ay: 


._ = . 29 Sin SHA Si 29 j = geee 5 i 1 
C; A. cos ong j j in oan ght 24] M, since the above equations 


then become 
(C, cos ano st, + S; sin ano st.) ent 
(6) 


M 
pa ano 5(C, sin ano st . S; cos ano st, ) a0 


which are linear in the new unknowns C\. C5. Sj. The A. and b can be 


recovered later from the C., 33 by means of the formulas 


au 
Ms 
capa wort. all 
J A J 
bey (7) 
2nd. = tan “J 
J C; 


We now review the reasons for prefering the Modified Gram-Schmidt least 
squares fit algorithm to that used in the harmonic analysis of hourly heights 
(Foreman, 1977). Given an overdetermined system of equations written in the 
matrix form Ax = b, a common approach to the linear least squares fit problem 
is to form and solve the normal equations AT‘Ax = A™b. These are not over- 
determined, but are frequently i11-conditioned, making the solutions very 
sensitive to roundoff errors, etc. In order to preserve accuracy, it is 
preferable to compute the least squares solution directly from the original 
overdetermined equations by orthogonalization procedures, such as Householder 
triangularization, singular value decomposition or the Modified Gram-Schmidt 
method. For instance, solution by normal equations requires double precision 
arithmetic to give the same accuracy as Householder's method achieves in 
Single precision (Barrodale et al., 1978). 


Nevertheless, when the number of equations N is much larger than the 
number of parameters M, the normal equations approach has its advantages. 
Not only do the formation and solution of these equations require about half 
as many operations as an orthogonalization technique, but if these equations 
can be formed directly from the data rather than from the overdetermined 
system, then only M(M+1)/2 storage locations, as opposed to at least NM, 
are required. This is, in fact, why the normal equations were used in the 
harmonic tidal heights analysis of hourly observations (Foreman, 1977). 
There, the normal equations are formed directly and efficiently (the use of 
trigonometric identities avoids rounding errors usually encountered with 
cumulative sums) and aggravation of the problem's already il1-conditioned 
nature is avoided. The storage savings can be significant - for instance, 
in a one-year analysis, where the approximate values for N and M are 8760 
and 137, respectively. In fact, on some installations, there may not be 
sufficient storage for an overdetermined array of these dimensions. 


Since the numbers of observations and constituents will generally be 
smaller when high-low analysis is used, storage considerations will be less 
important. Also, the identities used to form the normal equations in the 
regularly sampled case no Jonger apply. Orthogonalization methods are, 
therefore, preferable and the Modified Gram-Schmidt algorithm (Barrodale et al., 
1974) was selected, since it is competitive in all respects with the House- 
holder method and was already available on the Institute of Ocean Sciences' 
computer. 


Orthogonalization methods obtain the least squares solution to the 
matrix equation Ax = b by forming an equivalent system of equations which 
is easier to solve. In particular, the classical Gram-Schmidt technique 
does this by calculating an orthogonal set* of vectors 


(Sis elena such that Tor’k = Ty... fis the set {4,++-4, } spans the same 
k-dimensional subspace as the given set of linearly independent vectors 

{ay +++ak}, where the set {a,.--an+1} are the columns of the augmented matrix 
Atb arising from the overdetermined system Ax = b. The set of mathematical 
formulas which calculate the qs vectors iteratively are as follows: 


a etal (8) 
j-1 
ee Per ae 44 14 Pees Ki (9) 
F. 
a as qs 
where ie J i (10) 
q q 


+ the set of vectors {912+ +n} is said to be orthogonal if a4 °9, = 0 for 


ait i& 4. If, in addition, a; 4, = 1, the set is orthonormal. 


1] 


In order to convert these equations to matrix notation, let A be the matrix 
with columns ays Q be the matrix with columns q5 and R be the upper triangular 
matrix with unit diagonal elements and super-diagonal elements given by (10). 
Then equations (8) through (10) can be written as A = QR. 


The Modified Gram-Schmidt method is a variation of the classical 
technique which makes use of the fact that the value of the inner product 
asa; in equation (10) will not change if a, is replaced by any vector of 
the form 


a ey a oy Gy Cia 


where the a, are a set of arbitrary numbers. In particular, if one chooses 
the numbers a, SO as to minimize the norm of the vector ane it can be 
shown that replacing aj by aye in the classical Gram-Schmidt orthogonal- 
ization produces a method that is more stable numerically. This algorithm, 
which has become known as Modified Gram-Schmidt (MGS), is described by the 
following equations (Lawson and Hanson, 1974) : 


: reid 5 ee BT (12) 
Lael aj) (13) 
dinavFide 4 3} (14) 
3 ‘ aj!) a4 UR aligns oy El His 

ds? Po A Ue 
oF) 2 a). ng, | (16) 


In order to use the Modified Gram-Schmidt orthogonalization to minimize 
the sum of the squares of the residuals (i.e. || Ax - bil”) for the over- 
determined system Ax = b, first form the augmented m x (n + 1) matrix 
A' =[A:b]. A' is then orthogonalized to obtain 


Age =anQ' R' aa 
where the column vectors given by (13) constitute the m x (n + 1) matrix 


Q', and R' is an upper triangular matrix with unit diagonal elements and super- 
diagonal elements given by (15). Defining D' to be the (n + 1) x (Mead) 


diagonal matrix with elements specified by (14), a newm xm orthogonalt 
matrix Qo is introduced such that 


Ue (i = Q', where 0 has dimension (m-n+1) x (n+1) 


2 : 
This means that the first n+l columns of Qo will be q./d. fO%: th Ha) 1 
and the remaining m-n+1 need only complete the orthonormal set. 
Partitioning R' and D' into i Cc and |D 0 respectively, 
oe] 0 fen 


L 

where both D and R are n x n, we can.then write 
DR De 

Q 

re) 0 dn+1 


O #0 


Making use of the property that || B'x| lx | for any orthogonal matrix B, 


it then follows that 


Jax - bf] = a," (Ax - bd)” 
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v2 
JO(Rx-c) J + diay 


Therefore, the minimum value of IlAx - bll* is dua sf and is attained by the 
vector Xo which satisfied Rx = c. 


A later version of MGS (Barrodale et al., 1974)-is more 
efficient, in that the matrix equation Rx = c is solved as part of the 
orthogonalization process. Specifically, A' is now defined as the larger 


partitioned matrix A' = Ab » where I is the n x n identity matrix and 0 
£0 
is the n x 1 vector of zeros. Applying MGS to A' results in the following 
matrix of orthogonal columns | Q 0 where [Q r] is the Q' of equation 
Riper ee 
(17) and R,c are the same as in the partition of R'. 


Thus, the least squares solution Meee R-lc can be easily removed from this 


matrix. 
Peseta Pied fas SR Re Bi atte Mice ole Sdn, Dai 


+ An orthogonal matrix Q satisfies the condition q! =Q°?. AS’ a result, its 
columns form an orthonormal set. 


Moreover, since 


A b Q r Ria iG QR Qctr 
BGohen WO) cb et HRES) [Owed | vobk a) 0 


it follows-that, 0 = Naan b = ART ctr and hence r = b - AX, is the vector 
of residuals corresponding to the least squares solution of the over- 
determined system. 


5. MODIFICATIONS TO RESULTS OF LEAST SQUARES ANALYSIS. 
5.1 Nodal Modulation. 


The tidal potential contains many more sinusoidal constituents than 
are commonly sought in tidal analysis. Due to the small frequency differences 
between many of these (and the great length of record required for their 
separation) and the relatively small expected amplitude of some, it is not 
feasible to analyse for all of them. The standard approach is to lump 
together constituents which have the same first three Doodson numbers (see 
Godin, 1972) and to assume that each such cluster can be replaced by a 
single sinusoid having the same frequency as the major constituent (in 
terms of tidal potential amplitude) in the cluster. This major contributor 
lends its name to the cluster and lesser constituents are termed its 
satellites. An amplitude and phase is then calculated from the data for each 
apparent major constituent (in fact, for the replacement sinusoid representing 
each cluster). However, since these values represent the cumulative effect 
of all constituents in the cluster and the latter all differ slightly in 
frequency, the amplitude and phase of the replacement sinusoid vary Slowly in 
time and do not provide a basis for predicting the contribution of the 
cluster to the tidal signal at a subsequent time. To avoid this difficulty, 
the time-invariant amplitude and phase of the major constituent in each 
cluster are calculated from those of the replacement sinusoid. This 
adjustment procedure is known as “nodal modulation". To predict the 
contribution of the cluster at a later time, the major constituent is first 
calculated and then the nodal modulation corrections to its amplitude and 
phase are applied in the reverse sense to obtain the contribution of the 
cluster as a whole. 


t 


In this report, the replacement sinusoid for the j i cluster is first 


written as 


A; cos anfo;t, - 4; ] (18) 


where A. and ¢. are termed the raw amplitude and raw local phase. Time ty 
is measured from the mid-point of the record being analysed. A customary 
notation for (18) is 


fect) ka. 


Me os YT. 
: 5 608 ot [o5ty F Vid ges Me (19) 


J 


14 


which expresses the relation between the cluster contribution in terms of 

the amplitude, a; and phase, 9., of its major constituent. The nodal 
modulation terms of f(t) and u(t) vary slowly with time and, for records up 
to one year in length, very little error is introduced by assuming them to be 
constant and equal to their value at ty = 0, the midpoint of the record. 
Further details of the nodal modulation calculations and the evaluation of 
f(t) and u(t) are given in Foreman (1977, p. 42). 


5.2 Astronomical Argument Corrections. 


The astronomical argument correction arises from the need to express al] 
constituent phase lags with respect to a universal time and space origin. 
Instead of regarding each tidal constituent as the result of a particular 
component in the tidal potential, an artificial causal agent can be 
attributed to each constituent in the form of a fictitious star which travels 
around the equator with angular speed equal to that of its corresponding 
constituent. Making use of this conceptual aid, the astronomical argument of 
a given tidal constituent can be viewed as the angular position (longitude) 
of its fictitious star. For historical reasons, all such arguments or 
longitudes are expressed relative to the Greenwich meridian and can, 
consequently, be expressed as functions of time only. The replacement 
sinusoid (19) is often written as 


f(t) a cos 2r[ \; (t) + u(t) - 9, ] (20) 


where V.(t) is the longitude of the fictitious star relative to Greenwich 
and the new time variable t has an absolute datum such as some calendrical 


aio 


landmark. 9 is termed the 'Greenwich phase lag' of the j° constituent. 


Ignoring the minute effects of long-term changes in the astronomical 
variables (see Foreman, 1977, p. 15). 


VaCE, oF ° + ’ 
.{ ) a,t YO 


where V4 iS a phase correction due to the change of time datum. If te is 
the midpoint of the record on the new time scale, t, then t = t, 7 t and 


A, cos erlo st, - >: ] 


j a cos enfo.t wot = 6] 


J J 


= A. COS én[ V;(t) os LP i o jt. a ) 


= A; cos anf V(t) - V5 (te) - oJ 


fe 


Comparing this with (20), we see that the Greenwich phase lag is, in fact 


95 E >; ¥ V5(t.) + us(t.) 


5.3 Inference. 


Inference is the term used in tidal analysis to describe the extraction 
of certain important constituents excluded at the least squares fit stage on 
the grounds of insufficient record length but deduced afterwards from included 
constituents to which they bear a known amplitude and phase relationship. 
When accurate inference constants (amplitude ratio and phase difference) 
are available, inference not only yields amplitudes and phases for the 
inferred constituents, but also significantly reduces periodic variations in 
the estimated amplitudes and phases of the reference constituents. The 
computational steps involved in inference are given in detail in Foreman 
(1977). That material will not be repeated here, but should be read in the 
light of the following comments. 


The question of when constituents should be included directly in the 
least squares analysis, and when they should be inferred is not easily 
answered. The Rayleigh criterion, which is used to select constituents in 
the harmonic analysis of hourly tidal heights (see Foreman, 1977, aly Cog We Be 
incomplete in its presumption that a record of length T is required to 
distinguish constituents with a frequency separation of Te nt ract tt 
conflicts with the algebraic viewpoint whereby, for any four independent 
observations, one can obtain four equations and solve for four unknowns 
(2 amplitudes and 2 phases), regardless of the frequency separation. The 
missing consideration in both of these viewpoints is that sea-level 
observations contain, in addition to discrete tidal signals, contributions 
from a continuous noise spectrum of geophysical origin and from random errors 
in recording the observations. Taking these effects into account, Munk and 
Hasselmann (1964) showed that meaningful information can be gained about the 
frequencies o, and o» provided that 


| 0, —- 0} | > no 


ks taal Pf TRS 
(signal/noise level)* . 


It is interesting that essentially the same result can be derived by 
considering the sensitivity of solutions of a linear system to the condition 
number of its coefficient matrix. The following account is based on the 
detailed discussion in Ortega (1972). If K(A) = HANIA, is the condition 
number+ of matrix A and x,X are such that Ax = b and AX = b + 6b, then 


A 
Ix - xil 
Ix Il 


sbll (21) 

K(A) bt 

eee a a F Sap Pine BOO © nt A toot 9) te) ee Yel htt hoe 

+ The conventional condition number K(A) defined here differs from the norm- 
alized condition number C(A) calculated during the solution of the normal 
equations in hourly heights analysis (Foreman, 1977, p.38). Whereas K(A) is 
unity for a diagonal matrix and assumes higher values for more i11- 
conditioned matrices, C(A) lies between 0 and 1; 0 corresponds to a singular 
matrix and 1 to a diagonal matrix. 


In order to apply this result to the present problem, assume that 
Ax = b are the normal equations in matrix form arising from a least squares 
fit for the amplitudes and phases of several tidal constituents. In 
particular, assume that the right-hand sides b and b + Sb are respectively 
calculated from observations without and with background noise; that is, 
b assumes observations from a signal that is comprised purely of tidal 
components whereas b + 6b assume the same signal plus noise. The effect 
of seeking amplitudes and phases corresponding to frequencies o, and o» that 
are relatively close, i.e. Jog - oi] < T !, is to make the appropriate rows 
in A more linearly dependent (om Foreman, 1977, p.33, regarding the 
structure of A), and so increase K(A). Hence, in the presence of substantial 
background noise, one can expect a significant difference between the 
calculated set of parameters x and their “true” values, x. Assuming that 
accurate inference parameters are available, inference in such a case would 
yield better results, because solving for the parameters of only one 
frequency, o; OY oo, would remove the linearly dependent rows and reduce 
K(A). On the other hand, if the noise level is very small, the effect of a 
large condition number resulting from two close frequencies would be counter- 
acted and a reasonably accurate set of parameters x could be expected 
without inference. 


Table 2 gives the results of tests designed to demonstrate these points. 
Two 15-day records of hourly tidal heights were simulated, one using only 
the constituents Z,, 0;, Kz, Mz and So, and the other with these same 
constituents plus random background noise. Specifically, the tidal signal 
varied over the range [2.77, 47.23], while the uniformly distributed random 
noise added the range [-2.5, 2.5]. Three sets of six consecutive 60-hour 
harmonic analyses were executed; the first searching directly for all 
constituents; the second searching for only Z_, K;, and Mo; and the third 
one extending the second by inferring 0, and 4 from K, and Mo respectively 
(inference parameters were calculated from a 15-day analysis which sought all 
constituents). In order to compare performances ,means and standard 
deviations were calculated for each amplitude and phase over the six 
analyses in each series. 


Results from the analyses of the tidal record with no background noise 
(Tests 1 to 4, Table 2) demonstrate a clear advantage to seeking all 
constituents directly in the least squares fit. The small non-zero 
Standard deviations are attributable wholly to the fact that the data was 
rounded tofour digits, making Ilcbll slightly larger than zero. The standard 
deviations for the runs with inference (Test 4) are not zero because of 
simplifying assumptions in the inference method itself. 


However, when random noise with range [-2.5, 2.5] is also present in the 
tidal record (Tests 5 to 8), the standard deviations for the inference runs 
(Test 8) are consistently less than those obtained by the direct inclusion 
of all constituents in the least squares fit. This is a consequence of a 
reduction in K(A) from 120.1 for Test 6 to 2.884 for Test 7.+ 


+ The L. norm was used in equation (21). 
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(Corresponding values for C(A), the normalized condition numbers routinely 
included in the harmonic analysis program output, were 0.9697 and 0.0069. ) 
The average value of sbi for the six runs of Test 6 was 0.0531 while, 

for Test 7, it was 0.0479. Consequently, applying equation (21), we 
anticipate a max imum change in C. and $3 (s4) of 638% in the direct inclusion 
case and 13.8% when inference is used (assuming accurate inference constants 
are available). 


However, for most tidal analyses, equation (21) cannot be applied, since 
lob is unknown. Munk and Hasselmann (1964) derive a more useful formula 
for estimating the amplitude variances of close constituents in the presence 
of noise. Specifically, if the underlying noise spectrum is S(o) and the two 
neighbouring constituents have frequencies o, and oy, then the estimated 
variance of either amplitude is 3 S(o)n 2 | o2 - oO] ar Applying this 
result to the previous test data yields expected standard deviations of 0.576 
and 0.622 in the 0,/K, and M,/S» amplitude ratios, respectively. Although these 
values are much closer to those in Table 2 than the upper bound estimate 
derived from equation (21), it should be noted that they are all underestimates. 


For actual sea-level observations, Munk and Bullard (1963) estimate the 
geophysical noise level (due mainly to atmospheric excitation) at tidal 
frequencies to be S(c) = 1 cm2/(cycle per day). Applying this value to a 
30-day harmonic analysis in which both P, and K, are sought directly (these 
constituents require 183 days to differ by one cycle) yields an expected 
standard deviation of 0.613 cm in either amplitude. However, a sequence of 
twelve-monthly analyses at Prince Rupert, in which both P, and K, (and three 
other constituent pairs that also require approximately six months to differ 
by one cycle) was sought directly, produced standard deviation estimates of 
17.2 and 18.4 cm, respectively, for the amplitudes. Thus, in this case, 

S(o) = 1 cm2/(cycle per day) is a gross underestimate. It is worth noting 
that, with accurate inference parameters, the two amplitude standard 
deviations, after inference, reduce to 0.54 and 1.72 cm, respectively. (In 
each case, this is 3.3% of the amplitude.) 


PART III. COMPUTER PROGRAM FOR HIGH-LOW ANALYSIS. 
6. GENERAL DESCRIPTION OF PROGRAM, 


This program analyses irregularly-sampled tidal heights observations over 
a specified period of time. Although these observations are normally taken 
at high and low water, the program can also be used when the observations are 
not extreme values. Amplitudes and Greenwich phase lags are calculated for 
all requested constituents by a least squares fit method (84), coupled with 
nodal modulation (85.1). If the record length is such that certain 
important constituents cannot be resolved satisfactorily by including them 
directly in the least squares fit, provision is made for inference of their 
amplitudes and phases (85.3). 


6.1 Routines Required. 


Gl) MeN, <cze Gc reads input data, controls all output and calls 
other routines. 


(2 pees inna. & does a least squares fit (with the Modified Gram- 
Schmidt Algorithm) to find coefficients of the sine 
and cosine terms corresponding to each of the specified 
constituent frequencies. 


(Sovee VLEs ricrs ois reads required information and calculates the nodal and 
astronomical argument corrections for all constituents. 


fee acs chine reads required information and calculates the amplitude 
and phase of requested inferred constituents, as wel] 
as adjusting the amplitude and phase of the constituents 
used for inference. 


ChJ> COAY. . «date returns the consecutive day number from a specific 
origin for any given date in the range 1901 to 1999, 
and vice versa. 


Routines VUF, INFER and CDAY are identical to those of the same name in 
the tidal heights analysis program used for hourly observations 
(Foreman, 1977). 


6.2 Data Input. 


Three files or devices are used for input. File reference number 8 
contains the tidal constituent information that is necessary for the nodal 
and astronomical argument calculations; file inference number 9 contains 
the observed tidal heights and their times; and file reference number 5 
contains analysis type and tidal station information. A listing of file 
reference number 8, along with sample input corresponding to file reference 
numbers 9 and 5, are given in Appendices 3, 4 and 5, respectively. 
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File reference number 8 is a subset of the similarly-numbered file in 
Foreman (1977 §1.3). For most analyses, its contents should not require 
changing (see 6.4 for the circumstances under which this file might be 
changed). It contains the following three types of data, and is read in 
through entry point OPNVUF of subroutine VUF. 


(i) Two cards specifying values for the astronomical arguments 90.7 HO, 
PO, ENPO, PPO, DS, DH, DP, DNP, DPP in the format (5F13.10). 


SO -. mean longitude of the moon (cycles) at the reference time origin. 
HO - mean longitude of the sun (cycles) at the reference time origin. 


PO - mean longitude of the lunar perigee (cycles) at the reference 
time origin. 


ENPO - negative of the mean longitude of the ascending node (cycles) 
at the reference time origin. 


PPO - mean longitude of the solar perigee (perihelion) at the 
reference time origin. 


DS, DH, DP, DNP, DPP are their respective rates of change over a 
365-day period at the reference time origin. 


(ii) At least one card for all the main tidal constituents specifying their 
Doodson numbers and phase shift, along with as many cards as are 
necessary for the satellite constituents. The first card for each 
such constituent is in the format (6X, A5, 1X, 613, F5.2, 14) and 
contains the following information : 


KON = constituent name; 

II ,JJ,KK,LL,MM,NN = the six Doodson numbers for KON; 

SEMI = the phase correction for KON; 

NJ = the number of satellite constituents. 


A blank card terminates this data type. 


If NJ > O, information on the satellite constituents follows, three 
satellites per card, in the format (11X,3(313,F4.2,F7.4,1X,11,1X)). 
For each satellite the values read are : 


LDEL,MDEL, NDEL = the last three Doodson numbers of the main 
constituent subtracted from the last three 
Doodson numbers of the satellite constituent; 


PH = the phase correction of the satellite constituent 
relative to the phase of the main constituent; 


EE = the amplitude ratio of the satellite tidal potential 
to that of the main constituent. 


2] 


IR 1 if the amplitude ratio has to be multiplied by 
the latitude correction factor for diurnal 


constituents, 


= 2 if the amplitude ratio has to be multiplied by 
the latitude correction factor for semi-diurnal 
constituents, 


= otherwise if no correction is required to the 
amplitude ratio. 


(iii) One card specifying each of the shallow water constituents and the 
main constituents from which they are derived. The format is 
(6X,A5,11,2X,4(F5.2,A5,5X)) and the respective values read are: 


KON = the name of the shallow water constituent; 

NJ = the number of main constituents from which it is 
derived; 

COEF ,KONCO = combination number and name of these main 


constituents. 


The end of these shallow water constituents is denoted by a blank 
card. 


II. File reference number 9 contains only one type of data, the observed 
tidal heights and times. For convenience, the input format was chosen 
to be the same as the daily high-low output produced by the tidal heights 
prediction program (Foreman, 1977, p.54). Specifically, each record has 
the format (2X,15,213,12,6(13,12,F5.1)) and contains the following 


information: 

ISTN = tidal station number; 

LDZtM TY = day, month, and year of subsequent observations; 
ITH, ITM,HT = times (in hours and minutes), and height, of up to 


six observations for the specified date. If there 
are less than six observations for a day, they are 
padded to that number with the values 99,99 and 99.9 
for the times and heights, respectively. If there 
are more than six observations on a given day, as 
many records are included as necessary, each with 
the same repeated date. 


Missing days and/or missing observations (highs,lows) are permissible. 


However, it is necessary that the records be ordered according to date. 


LIL, 


La 


File reference number 5 contains five types of data. 


(i) 


(ii) 


Ltn 


(iv) 


One record for the variables MF,IDERV,WT in the format (215,F5.2). 


MF = the number of constituents, including the constant 
term, Z,, to be included in the least squares fit. 


IDERV = ] if all observations are.extreme values and it is 
desired to use the derivative conditions in the 
least squares fit, 

= 0 otherwise. 

WT = the weight to be applied to the derivative 
condition when IDERV = 1. A recommended value is 
Wh Suke0. 


One record for each of the MF constituents to be included in the 
fit. Each record contains the variables NAME and FREQ in the 
format (A5,2X,F13.10). NAME is the constituent name, which should 
be left-justified in the alphanumeric field, while FREQ is its 
frequency measured in cycles per hour. In order that there be 
sufficient information available to calculate the astronomical 
argument and nodal corrections, all these constituents must be 
included in the list given in Appendix 2. The order in which the 
constituents are input is also the order in which the results are 
output. The constant term Z, must be first. 


One record in the format (6(3X,12)) containing the following 
information on the time period of the analysis: 


ID] ,IMI,IY1 - day, month, year of the beginning of the analysis 
period; 
D2, iM2.1y2 - day, month, year of the end of the analysis period. 


One record in the format (15,5A4,1X,A4,415) containing the 
following tidal station information: 


JSTN = tidal station number; 
(NSTN(I),I=1,5)= tidal station name; 
LTA 


i} 


time zone in which the observations were recorded; 
LATD,LATM = station latitude in degrees and minutes; 


LOND ,LONM = station longitude in degrees and minutes. 
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(v) One record for each possible inference pair. The format is 


(2(4X,A5,£16.10) ,2F10.3) and the respective values read (through 
entry point OPNINF of subroutine INFER) are: 


KONAN and SIGAN = the name and frequency of the analysed 
constituent to be used for the inference; 


KONIN and SIGIN = the name and frequency of the inferred constituent; 
R = amplitude ratio of KONIN to KONAN;: 


ZETA = the Greenwich phase lag of the inferred 
constituent subtracted from the Greenwich phase 
lag of the analysed constituent. 


These are terminated by one blank record. 


As before, constituent names should be left-justified in the 
alphanumeric field, frequencies are measured in cycles per hour, 
and all constituents must belong to the list in Appendix 2. 


6.3 Output. 


At present, only line printer (file reference number 6) output is 
produced by the high-low analysis program. The first page simply echoes the 
requested constituents to be included in the fit, the analysis period and 
tidal station information, and the heights and times of the observed heights. 
The second page notes whether or not the derivative conditions were used for 
all observations and, if so, the value of the requested weight. It also gives 
the following direct results from the least squares fit: Z, amplitude and 
coefficients of the cosine and sine terms of all other constituents; the 
largest residual value in the overdetermined system of the equations and the 
residual sum of squares; the standard deviation of the right-hand sides of 
the overdetermined system and the root mean square residual error. The third 
and final page gives the raw amplitudes and raw local phases followed by the 
nodally-corrected amplitudes and Greenwich phases for all constituents. If 
there has been inference, these values are repeated with correction, and the 
new residual root mean square value is specified. 


In the event that a least squares solution cannot be found because a 
dependent column is encountered during the orthoganalization procedure, a 
message to this effect along with the suggested corrective procedure is 
printed. If the column dependency is borderline, a slight increase in the 
value of the variable TOLER may be sufficient to obtain the solution. 

However, it is generally better to remove from the least Squares analysis, the 
constituent (or its nearest neighbour, depending on which one has the smailer 
expected amplitude) corresponding to this tidal coefficient parameter. That 
is, if column 2n-1 or 2n is dependent, then remove constituent n or its 
nearest neighbour. If inference parameters are available, the amplitude and 
phase for this constituent can still be obtained indirectly through inference. 
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The final page of output produced by the sample data input found in 


Appendices 3, 4 and 5 is listed in Appendix 6. 


6.4 


Program Conversion, Storage and Dimension Guidelines. 
The high-low analysis program was tested on the UNIVAC 1106 computer 


installation at the Institute of Ocean Sciences, Patricia Bay. Although the 
program was written in basic ASCII FORTRAN, some changes may be required 
before it can be used on other installations. These may include: 


(i) 


(iii) 


replacing all calls to routine INPROD in subroutine MGS when the FORTRAN 
compiler does not permit a single column of a two-dimensional array to 

be passed to a one-dimensional array through a subroutine call. Such 
changes will not be necessary when FORTRAN compilers store two-dimensional 
arrays by columns (and this is the standard FORTRAN convention). 

However, if this condition is not met, the INPROD calls are located in 
lines 102, 122, 135, 156 and 173 of subroutine MGS and the replacement 
code is specified in the comment statements preceding these CALL 
statements. 


deleting some or all REAL*8 declaration statements. These are used 
either to permit alphanumeric strings longer than 4 characters (which is 
the ASCII FORTRAN single precision word-length) to be used (e.g. 
constituent names are read in A5 format), or to gain computational 
precision in the average, standard deviation, and vector inner product 
calculations. CDC installations may not require double precision 
variables for the latter of these cases because their single precision 
accuracy is about 1.5 times that for IBM and UNIVAC (= 8.1 digits), and 
do not require them for the former because their single precision word- 
length of six characters is longer than the maximum alphanumeric string 
Size of 5 required by this program. 


altering the variable list structure for the ENTRY statements OPNINF 
and OPNVUF, and references to them. 


changing some, or all of the file reference (or device) numbers from 
their present values, in order to conform with local machine restrictions. 


altering the input tolerance, variable TOLER, for the Modified Gram-Schmidt 
routine. If the inner product of an orthogonalized column with itself 

is less than TOLER, the column is considered to be dependent. Typically, 
TOLER is chosen to be less than 10**(-D) where D represents the number 

of decimal digits of accuracy available. However, if the overdetermined 
matrix 1s poorly scaled, it may be necessary to either choose a much 

larger value or remove the corresponding constituent from the analysis. 

On the UNIVAC, a conservative value of 10**(-7) was chosen for TOLER. 


The program, in its present form,requires approximately 3000 and 73000 


UNIVAC words for the storage of its instructions and arrays, respectively. 
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A large part of this is for the array Q which stores the overdetermined 
system of linear equations and is presently dimensioned to handle approx- 
imately 800 observations with the derivative condition and 20 constituents. 
If the analysis is much smaller than this and memory requirements are 
restrictive on a particular installation, or there is a need to economize, 
the program size can be cut significantly by reducing the size of this array 
and resetting variables NMAXP1 and NMAXPM appropriately. 


In the event that changes are required to the program, restrictions on 
the minimum dimension of all arrays and minimal values of special parameters 
are as follows: 


Let 


MC be the total number of constituents, including Z, and any inferred 
constituents; 


NOBS be the number of tidal height observations; 
NR be the number of input records of observed tidal heights. 
MPAR be 2*MC - 1. 


NEQ be NOBS*2 if all the observations are extremes and the derivative 
condition is to be included for each, and NOBS otherwise. 


Then, in the main program, parameters MXNDAY, NMAXP1 and NMAXPM should 
be at least NR, MPAR + 1, and NEQ + MPAR respectively; arrays FREQ, NAME, 
AMP, PH, AMPC, and PHG should have minimum dimension MC; arrays X and Y 
should have minimum dimension NOBS; arrays ITH, ITM, and HT should have 
minimum dimension 6; array P should have minimum dimension MPAR; and 
array Q should have the exact dimension of NMAXPM by NMAXP1. In sub- 
routine MGS, arrays Q and X have variable dimensions, and should be the 
same as Q and P in the main program. 


The dimensions of any array in subroutine VUF need only be changed if 
a new constituent is to be added to the list in Appendix 2. In such an 
event, the contents of the data file associated with file reference number 
8 must also be augmented in order to permit calculation of astronomical 
argument and nodal corrections for this constituent. In order to under- 
stand the structure of this file and the resultant calculations, consult 
Foreman (1977). Restrictions on the minimal array dimensions can be found 
there, also, aS well as in the comment statements of the subroutine itself. 


In subroutine INFER, array KON is passed in the argument list from 
the main program and so need only be dimensioned 2; and arrays KONAN, 
SIGAN, KONIN, SIGIN, R and ZETA can presently accommodate a maximum of 9 
inferred constituents. 


In subroutine CDAY, arrays NDP and NDM should have dimension 12. 
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Abstract 


Oceanographic observations were taken through the ice in Robeson 
Channel during April and May 1971. These observations were water bottle 
samples for temperature and salinity, current meter profiles to depth of 
70 metres and continuous current meter observations under the ice from 
each side of the channel. | 
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During the spring months of April and May 1971, oceanographic 
observations were made through the ice in Robeson Channel at locations 
Situated between Ellesmere Island and Greenland (Figure 1). These 
observations were part of a cooperative program carried out by Defence 
Research Establishment Pacific of the Department of National Defence and 
Marine Sciences Directorate of the Department of Environment (now Ocean and 
Aquatic Sciences, Department of Fisheries and Oceans). 


PROGRAM: 


First, a main-camp tent was set up on the ice in mid-channel of 
Robeson Channel opposite Wrangel Bay (Figure 2). A hole was then blasted in 
the 1 metre thick sea ice and a work tent, complete with a gas-fired heater, 
was erected over the hole (Figures 2 and 3). A light bulb, powered by the 
main camp electric generator (when operative) was inserted into the ice hole 
to prevent it from re-freezing, however, daily chiselling was also required 
to loosen newly formed ice so it could be flushed out of the hole. 


At the main camp, three oceanographic casts and several bathythermo- 
graph lowerings and plankton hauls were made through this hole. A continuous 
recording current meter was suspended on a wire on one side of the hole but 
was later found to have malfunctioned. A time-series of vertical current 
velocity profiles below the ice was taken with only some success due to 
malfunctions of the current meter's direction sensing unit. With the use of 
a helicopter, continuous recording current meters were transported to two 
additional locations on both sides of Robeson Channel and installed at 
5 metres' depth below the ice (Figure 4). These were in operation for 15 
days. Two oceanographic casts were also made at these two current-meter 
Sites. 


Equipment Used 


Recording Current Meters. A Braincon current meter with a photo- 
graphic film recorder was used with a modified collapsible tail (Figure 4). 
This tail folded in against the body so that the assembly could be lowered 
down a 9" hole. A candy "life-saver" held the tail into the body of the 
meter and after approximately 20 minutes the "candy" dissolved and released 
the tail to its operating position under the ice. When the current meter was 
retrieved, the tail assembly was forced to fold downwards as the meter was 
pulled up the hole. 


Current Profilers. A Hydro Products Savonius rotor current meter 
was uSed to obtain a number of vertical current profiles. Some trouble was 
experienced with the direction sensing unit which depends on a spring-loaded 
contact on a toroid resistor. Also, the Savonius rotor jammed on several 
occasions due to a build-up of ice crystals on its reed switch. For longer 
periods of observation a Rustrak recorder was used during the silent hours 
and when no one was up. 


Hole Stoppers. "Weiner"-shaped rubber bladders (Figure 4) were 
inflated in the ice holes on both sides of the channel to prevent the holes 


from freezing in. The bladders leaked slowly from their valve stems causing 
the bags to soften, resulting in partially refrozen holes. 


Water Sampling Bottles. Fjarlie water bottles were used to collect 
water samples because of their ability to pass easily through a 9" diameter 
hole. The bottles were housed in a heated, insulated, 6-bottle carrying case 
for use from the helicopter. 


Sample Bottles. Samples were drawn into 5-ounce citrate bottles. 


Bathythermographs. A bathythermograph, operable to a depth of 
270 metres, was used and functioned well provided that the unit was warm prior 
to inserting it into cold water. 


Oceanographic winch. A winch, designed and built at the DREP 
workshop, was used for raising and lowering oceanographic instruments. This 
winch (Figure 3) worked well when driven by a 3/4 horsepower electric drill 
motor run off the main a.c. power generator. The unique design of the winch 
and block assembly, which stood above the hole, permitted one man to operate 
it from inside the tent. When the winch was used away from the main camp 
and transported by the helicopter, a Stihl power-saw motor was used to 
provide power. 


OCEANOGRAPHIC OBSERVATIONS 


Each of the five oceanographic casts, three at the main camp and 
one on each side of Robeson Channel, sampled temperatures and salinities vs 
depth. The samples were obtained using reversing water bottles and 
thermometers. After the reversing bottles were retrieved and when the 
auxiliary thermometer had reached equilibrium, the reversing thermometers were 
read and recorded. The samples for salinity determination were drawn into 
5-oz citrate bottles and stored for later analysis at the Pacific Environment 
Institute using an inductive salinometer. 


Figures 5 to 9 inclusive are plots of salinity (S), temperature 
(T) and sigma-T (G) vs depth in metres. Locations are shown in inset in 
Figure 1. STN-ROB1, STN-ROB2 and STN-ROB3 were at the main camp in the middle 
of Robeson Channel, STN-ROB4 was on the Greenland side, and STN-ROB5 was on 
the Ellesmere Island side. (The data from which these plots in Figures 5 to 
9 inclusive were obtained is listed in Appendix I.) Figure 10 compares cross- 
channel salinities and temperatures in the spring of 1971 (upper two diagrams) 
with data collected in the same section in the summer of 1971 (Herlinveaux and 
Sadler, 1973) (lower two diagrams). Figure 1] shows temperature vs salinity 
(TS) plots for all five oceanographic casts taken during the spring of 1971. 


The bathythermograms are presented in Appendix II and are hand- 
drawn reproductions of the slides. 


Currents 
Current meter records, each of 15 days' duration and taken from a 


depth of 5 m at each side of Robeson channel (Figure 1), were digitized and 
processed on a computer. (The raw data print-out for the two positions is 


listed in Appendix III.) Frequency distributions of direction and speed are 
shown in Figure 12 for the east side of the channel and Figure 13 for the 
west side. Polar coordinate histogram plots of direction and direction 
frequency are shown in Figures 14 and 15. Vector trajectory plots of water 
movements on either side of the channel are plotted in Figures 16 and 17. 
Figure 18 shows time-series plots of speed and direction at both locations. 


The current velocities obtained at the mid-channel camp at 
various depths down to 100 metres are presented, in part, in Figure 19. 
There were occasions when current meter directions were erratic and therefore 
unreliable. The currents shown in Figure 19 were selected to give the reader 
an idea of the variation with depth observed. The rest of the data remain 
on file. 


When not being used for profiling, the current meter at the camp 
was fitted with a two-channel Rustrak recorder and suspended at various 
depths. These data are not presented but are on file at the Institute of 
Ocean Sciences, Patricia Bay. 


INTERPRETATION 


The oceanographic data shown in Figures 5 to 9 inclusive show a 
50 metre thick upper layer above a halocline and thermocline extending down 
to 200 metres. Below 200 metres, there is a gradual increase in salinity 
and temperature. Minor temperature inversions occur in the upper 30 metres 
probably caused by cloud-like water layers moving in the channel. Most of 
the variation in the water column occurs at and below sill-depth of 130 
metres, which limits the exchange of deeper-water into the Nares Strait 
system. A comparison of Figures 5, 6 and 7 will give some idea of the 
variation to be expected over a period of 13 days. 


The cross-section distributions of temperature and salinity 
observed in the spring and summer of 1971 (Figure 10) show that isotherms and 
isohalines tend to slope upwards to Ellesmere Island in the upper 100 metres 
and downward below 200 metres. In the upper layer, temperatures from spring 
to summer increased by about 0.7°C at the surface and by about 0.2°C at 
50 metres. At mid-depth, very little change is seen and, at the bottom, the 
temperature decreased by 0.1°C. In summer salinities decreased at the 
surface possibly due to runoff, ice melt or both. In both spring and summer, 
the salinities were least along the Greenland side, probably associated with 
runoff from Greenland drainage. This suggests there is a component of flow 
to the north or at least a more sluggish southerly flow both in spring and 
summer. The slopesof the isotherms and isohalines further suggest that there 
are differences in the channel flow normal to its cross section - likely more 
vigorous in the down-channel direction on its westerly side. 


Figure 11 is a TS plot derived from the five oceanographic casts. 
The plots are arc-shaped with the centre of the curve at about 34.00%, which 
corresponds to the lower end of the halocline. This TS suggests there are 
only two water masses involved in this area - warm bottom water from the 
Atlantic, via Baffin Bay, and surface Arctic Ocean water moving south from 
the Lincoln Sea. 


Figures 14 and 15 are polar plots of all the current directions 
and speeds recorded from the west side and east side of Robeson Channel 
respectively. These figures show the number of observations in each 10° 
segment of direction and are an alternate way of presenting the data shown in 
Figures 12 and 13. 


To illustrate the net movements in Robeson Channel, vector trajec- 
tory plots are shown in Figures 16 and 17. Comparison of the two figures 
shows that, for the 14 days, the net movement down-channel was 17 km on the 
east side and 50 km on the west side. The down-channel residual currents are 
three times greater on the west side than the east side at 5 m depths - an 
expected result of the Coriolis effect acting on an otherwise uniform cross- 
channel flow. : 


Figure 18 illustrates the tidal influence on water movements. It 
is evident that there is a marked semi-diurnal fluctuation in current 
direction and amplitude, more so during the spring tides but tending to be 
diurnal during neap tides. 


The water movements in Robeson Channel have a strong tidal 
component - more in evidence during spring tides than during neap tides. 
Strong southerly trending residual currents are evident in Robeson Channel 
which, at times, exceed the magnitude of tidal currents. Therefore the best 
time to work in Robeson Channel on the bottom or on the foreshore would be 
during the neap tides. 


Figure 19 shows the current directions and speeds from the surface 
to 80 metres, measured in mid-channel. There were many profiles obtained, 
but the direction unit malfunctioned at times although the speed unit did 
work. Figure 19 shows that the currents generally increased with depth down 
to 80 metres. The swiftest currents observed were 50 cm/s at a depth of 
75 metres, while movement immediately under the sea ice was 15 cm/s. 


It is recommended, before any further physical oceanography is 
planned for Robeson Channel, that the data in this report should be analyzed 
in concert with that reported by Chow, 1975 and Sadler, 1976. 
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Speed and direction plots for the two sets of data collected on 
the east and west side of Robeson Channel. 


Examples of three speed and direction plots with depth observed 
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Figure ] Location of Robeson Channel south of Lincoln Sea and 
the positions of oceanographic observations. 


Figure 2 Ice Camp layout mid channel Robeson Channel with living quarters 
and electronics tent on left and oceanographic tent on right. 
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Figure 3 Inside view of oceanographic tent with specially designed 
winch. 
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Figure 4 Current meter being retrieved from ice hole with rubber tube 
hole stopper laying on ice. 
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Figure 6 Temperature, salinity and sigma-T distribution at 
mid-channel Robeson Channel. 
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Figure 7 = Temperature, salinity and sigma-T distribution at 
mid-channel Robeson Channel 
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Figure 8 Temperature, salinity and sigma-T distribution at 
the east side of Robeson Channel. 
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Figure 9 Temperature, salinity and sigma-T distribution at 
the west side of Robeson Channel. 
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Figure 19 Examples of three speed and direction plots with depth 
observed through the ice at the mid-channel camps. 
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ABSTRACT 


The copper and zinc concentrations in sediments of Georgia Strait, 'h..c., 
adjacent to the inactive iron-copper mine on Texada Island, are reported. 
Copper concentrations in 32 samples from 16 ug woe to 1584 ug gol and exhibited 
a distribution pattern relating these elevated concentrations to the mine as 
source. No significant elevation of zinc in the samples (range 28-195 ug gol, 
dry weight) was observed although a statistically significant relationship 
between copper and zinc levels was established. | 


INTRODUCTION 


Texada Mine, located on the northwest shore of Texada Island in Georgia 
Strait, B. C., closed down operations in September 1976. Since 1952, it had 
produced iron and copper concentrates for Japanese buyers. Until 1964, it 
was operated as an open pit. From 1966 on all operations were underground. A 
high grade of iron concentrate averaging 65% iron was produced. Copper concen- 
trates averaged 0.055 percent. Usable amounts of sulfur and phosphorus were 
also mined. 


Two factors were instrumental in bringing about the shut-down in 1976. 
The first and major one was the expiry of the Japanese contract on September 
30. Secondly, the mine would have been required to install a submarine dis- 
charge pipe to eliminate pollution of neighboring beaches. The existing dis- 
charge line consisted of a wooden flume which terminated above the high-water 
line. Periodically, the flume was extended as the tailings accumulated on the 
beach. 


Studies in 1971 (Goyette, personal communication, 1979) indicated that 
oysters (Crassostrea gigas) from the adjacent beaches were heavily contami- 
nated with copper. Concentrations averaged 1,200 ug g~! (range 1,128-1,240 
ug g~!) for the specimens. The condition of the few oysters found in this area 
was extremely poor. 


Discharge of solids to Georgia Strait averaged 3,266 tonnes per day (3,600 
short tons) in 1970. Because of the high limestone content, a good deal of the 
tailings were used for cement manufacture and local road construction. There 
exists at present a large limestone deposit on the property which the company 
is now utilizing for similar purposes. 


In January 1975, a cruise to collect samples of sediments on a grid in the 
Strait of Georgia was undertaken. This was part of a continuing program in- 
tended to study the chemistry of deposited tailings from various mines on the 
B. C. coast (Thompson and McComas, 1974; Thompson and Paton, 1975, 1976, 1978). 
We report here the data obtained for copper and zinc determinations for the 
32 sediment samples that were collected. 


METHODS 


Sampling 


Sampling stations and the area of study are shown in Figures 1 and 2. All 
samples were obtained with a Hydro Products Shipek grab. A vertical subsample 
was removed from the grab bucket with a plastic scoop, placed into a "Whirlpak" 
plastic bag and frozen immediately. Sampling was carried out from the C.S.S. 


Vector. 


Analysis 


In preparation for analysis the entire sub-sample was freeze-dried, 
pulverized in mortar and pestle and sieved through a 149 wm (100 mesh) screen. 


Analysis was made only on the < 149 ym fraction. 


An accurately weighed one gram sample was treated with hot aqua regia as 
described previously (Thompson and McComas, 1974). After dilution to 100 mL 
the extract was analyzed for copper and zinc by flame atomic absorption spectro- 
photometry. The instrument employed was a Jarrell-Ash, model 810. Instrument 
settings were primarily those recommended by the manufacturer. Data for copper 
and zinc were calculated as ug g7! of dry sediment. Analytical precision was 
determined by replicate (10 analyses of U. S. National Bureau of Standards 
copper mine tails [SRM 331]). Means of 903 + 6 (lo) ug.g~! and 36.3 + 5.9 (lo) 
ug.g 1 were obtained for copper and zinc respectively. The certified value 
tor Cu is 910 2.10 tis) ug.g7!. Precision for three replicate determinations 
of copper in Texada samples was 16.4 + 2.3 (lo) ug g-! (Station 28) and 1076 
+ 21.9 (1c) ug g-! (Station 7). 


RESULTS AND DISCUSSION 


Copper and zinc concentrations for the 32 sediment samples are shown in 
Table 1. Spatial distributions have been plotted as isopleths shown in Figures 
3 and 4% 


Copper concentrations are highest in the vicinity of the mine outfall, 
reaching a maximum of nearly 1600 ug g !, dry weight at Station 4. From the 
isopleth diagram it is apparent that there was considerable alongshore drift 
northwestward. These values are the highest for any mine-derived materials 
studied by us. Southeast alongshore drift is indicated also, however data 
in this direction are not sufficient to indicate the extent of this distribution. 


The structure of the isopleths in Figure 3shows that a second transport 
component lies directly south of the mine. This component probably accounts 
for a major proportion of the waste materials which were released in the 24- 
year period of operations. The outer edge of mine-derived material detectable 
by chemical means lies between the 50 and 100 yg g-! lines. A survey carried 
out by contractors for the mine in 1971 (Goyette, personal communication, 1979) 
delineated the distribution based upon visual observations. The outer limits 
as estimated have been reproduced in our scale in Figure 5. It can be seen 
that the direction of transport coincides well with our findings although as 
might well be expected the limits are one-quarter to one-fifth those determined 
chemically. 


Figure 4 isopleths indicate no discernible trends in zinc distribution, 
and in particular, there are no contours which suggest that the mine was a 
source of this element. The copper-zinc ratio has been used previously 
(Thompson and McComas, 1974; Thompson and Paton, 1975) as an operational para- 
meter to delineate mine tailings distribution. As an alternative to this 
graphic approach, which now appears to have limited applicability, we have 
plotted copper concentrations against the zinc values for each of the 32 sta- 
tions (Figure 6). We obtained two regression lines using reduced major axis 
calculations (Till, 1974) which eliminate biases imposed by normal least squares 
regression. They are used where there is no apparent dependence of one 
Parameter upon another. Inspection of the plot of all 32 data pairs indicated 
that stations 1 to 15 represented a different population that stations 16 to 32. 
This coincided roughly with the apparent zone of influence of the tailings 


(Figure 3). After establishing a normal distribution of data by x2 applica- 
tions two reduced major axis lines were generated (Figure 6) from data for 
stations 1-15 and 16-32. Both lines had highly significant correlation co- 
efficients of r = 0.57 (p < 0.02) and r = 0.90 (p << 0.001), suggesting a 

strong geochemical relationship between these two elements in both mine-derived 
and natural sediments (i.e. Stations 16-32). Although there is strong evidence 
to suggest that these two data sets do represent different populations there 
were not sufficient grounds to warrant a comparison of slopes which would serve 
as a statistical demonstration that deposits of material different from naturally 
occurring sediments exist in Georgia Strait adjacent to the mine. 


Most British Columbia coastal mines are located in fjords or inlets where 
water circulation is characteristically restricted. In cases where tailings 
have been released directly to the receiving waters (e.g. Howe Sound, Rupert 
Inlet, Alice Arm) the wastes have accumulated within a few km? of the discharge 
point. This type of situation in turn has created a potential or real source 
of heavy metal contamination due to diagenetic processes which serve to release 
the metals to interstitial and overlying waters. In the upper basin of Howe 
Sound, for example, we have observed enrichment of the water column and inter- 
stitial waters with copper (Thompson and Paton, 1978). There is some evidence 
also (Thompson, unpublished data) to suggest that pore waters in tailings 
deposited in Rupert Inlet may contain elevated concentrations of dissolved 
arsenic species because of the high arsenic content in the discharged materials. 


Texada Mine and the small, now defunct operation at Jordan River on 
Vancouver Island (Ellis, 1975) share the distinction of being on open, well- 
flushed coastlines. This situation has served to distribute waste materials 
over a wide area, which might conceivably produce two opposing effects. On 
one hand this dispersal would prevent, except in the vicinity of the outfall, 
an accumulation of several centimeters of tailings which might be required 
before discernible chemical impact. However, on the other hand, if contaminants 
contained in the tailings were particularly mobile this dispersal may have the 
effect of distributing the contaminant over a larger area. Population diversity 
and the concentrations of some heavy elements (including both metals and metal- 
loids) in organisms near industrial activity are often used as indicators of 
metal contamination. There are surprisingly few data available to indicate 
whether shoreline benthos on Texada Island have been affected. Goyette (1975) 
reported that oysters (Crassostrea gigas) collected within the tailings plume 
while the mine was still operating contained amounts of copper averaging 3,700 
mg kg-! (presumably dry wt.) and invertebrates dropped from 40 species in 
control areas to 2 species within the deposition area. 


There are insufficient data however to indicate the extent of contamination 
in intertidal or subtidal invertebrates. 


The entire question of metal concentrations derived from mining activities 
in benthos on this coast is one which begs an answer. It appears as if there 
may be a certain degree of site-specificity involved. This is suggested by 
data for oysters in particular. In Howe Sound this species (Goyette, 1975) 
contained up to 2,550 ug gol (dry wt.) south of the mine outfall, one year 
after operations ceased. Concentrations are considerably greater than those 
found for any pristine inlet. Although the evidence suggests strongly that 
the Britannia Mine activity was the source of this excess copper there is a 
possibility that the contamination is the result of leaching from natural 


mineral formations which contain the metal. 


Waldichuk recently (1979) discussed the differences in metal content of 
biota from two dissimilar mining activities. Nine years of monitoring have 
yet to show evidence of increased tissue burdens of copper or zinc in shellfish 
in the Rupert/Holberg Inlet areas. However there were considerable increases 
of zine in just one year of mining operation in mussel (M. Edulis) tissue from 
shore stations near the Marmorilik mine in Agfardlikavsa fjord in Greenland. 
Output of the latter mine is only 1400-1800 tonnes d~! compared to 38,000 
tonnes d-! for the Island copper mine. The nature of mineralization is 
different, which may allow for a greater solubilization of lead and zinc. 


Although there are no planned monitoring activities for the shoreline area 
of Texada Island it would nonetheless be of interest to see whether or not 
further uptake of copper by shellfish occurs in this area in spite of the non- 
restricted water movement. 


CONCLUSIONS 


The copper and zinc concentrations in marine sediments from Georgia Strait 
in the vicinity of a copper mine on Texada Island have been determined. Con- 
centrations of up to approximately 1600 ug g~! (dry wt.) were found in nearshore 
stations within 1.8 km of the outfall. Concentrations decreased rapidly with 
distance from the mine, with background values being observed at about 9-18 km 
distance. Zinc concentrations showed no trends with distance. 
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Figure l 
The coastline of southern British Columbia. Hatched area indicates 
the location of the study. 


Figure 2 
Study area showing station locations. 


Figure 3 
Isopleths indicating copper concentrations in sediments in ug g 
(dry weight). 


1 


Figure 4 
Isopleths indicating zinc concentrations in sediments in ug g. 
(dry weight). 
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Figure 5 
Delineation of tailing distribution obtained by visual identification, 
1971 (Goyette, private communication). 


Figure 6 
Reduced major axis regression lines for copper vs zinc in sediment 
samples. Equation for --(J-- Cu = 13.2 Zn - 937; r = 0.57. 


--@-- Cu = 1.00 Zn - 33.0; r = 0.90. 


Table 1 


Copper and Zinc Concentrations in Shipek Grab Samples 
of Georgia Strait Sediments near Texada Mine 


Station Cu Zn 
ug gw? ug go} 

1 82 32 
Z 88 28 
3 111 37 
4 1584 195 
5 1426 170 
6 370 73 
7 1076 130 
8 714 130 
9 265 90 
10 504 148 
Li 89 92 
12 281 106 
13 279 142 
14 187 104 
15 238 139 
16 139 146 
17 EZ 138 
18 105 143 
19 66 125 
20 85 138 
21 42 88 
22 95 139 
23 68 . 121 
24 32 43 
25 /M| 45 
26 DS 64 
2/ 44 76 
28 16 28 
29 asl, 89 
30 58 98 
31 81 59 
32 26 56 
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I. SUMMARY AND CONCLUSION 


This report summarizes the progress and results of the 1978-1979 studies 
contracted in support of Ocean Dumping legislation and reviews a workshop on 
ocean dumping research held on May 14, 1979 at the Institute of Ocean Sciences, 
Patricia Bay, Sidney, B.C. 


Alberni Inlet has been the focus for extensive ocean dumping studies 
since the registration in Canada of the Ocean Dumping Control Act, 1975. The 
Inlet serves as an ideal location for a case study in that it receives . 
relatively few other wastes, and has been studied intensively for its fisheries 
and for the environmental impact of a pulp mill and saw mills near its head. 
General explorations were made in the 1976-77 studies to obtain some 
understanding of the effect of disposal of wood wastes on the inner basin. 
The 1977-78 studies dealt with selected topics of dredge spoil impact on the 
Inlet with more emphasis on the oxygen budget of deep waters and the flux of 
oxygen to the sediments. 


The Ocean Dumping Control Act specifies allowable limits for certain 
contaminants in materials to be dumped. Several contracts were designed to 
establish meaningful allowable levels of contaminants specifically applicable 
to the Pacific Coast marine environment. The 1978-1979 fiscal year saw the 
continuation and completion of these contracts and as a result, there has 
been a general progression away from using Alberni Inlet as a focal point for 
ocean dumping research. 


Field studies investigating seasonal variation in the biological 
communities at the Port Alberni dumpsite have been completed. A preliminary 
review of the results demonstrates an interesting comparison between the 
control station and the dumpsite station. During the winter months and well 
into spring, when seasonal dumping occurs, the dumpsite contained fewer species 
and lower biomass than the control station. After May, and after dumping has 
ceased, the diversity at the dumpsite almost doubled; in some months the 
dumpsite had more species than the control station. Numbers of individuals 
at the dumpsite also increased in late spring and summer, but the biomass of 
the two communities remain relatively low throughout the summer period. It is 
apparent, that the dumpsite is quickly recolonized after dumping by such 
opportunistic species as Axinopsida serricata, Capitella capttata, Dorvillea 
spp., etc. Statistical analyses (contracted for 1979-80) will support a more 
definitive discussion of results. 


A natural history study has been done on various species predominating 
the four "biotic provinces" as separated statistically from the 1976 dumpsite 
survey. Feeding strategy and mode of movement are the most available data in 
the literature. Two species of polychaetes, Capitella capttata and 
Schistomeringos longtcornis inhabit the "dumpsite province" and exhibit life 
history characteristics which enable rapid colonization in degraded habitats. 
Wood-burrowing organisms, the mollusc Xylophaga washingtona, and the isopod 
Linmorta ltgnorun, were the prominent species in the "wood debris province". 


Physical parameters influencing the renewal of oxygen to the deeper 
waters of the inner basin comprised the remaining topics of research done in 
Alberni Inlet. Evidence available suggests that upwelling off the west coast 
of Vancouver Island is a response to the relaxation of south-east winds. 
Upwelling of sufficient strength tO induce renewal of water to the inner basin 
of Alberni Inlet is likely to occur every year. Consequently, deleterious 
effects due to oxygen demand should have no long-term accumulation. A 
mathematical model has been designed which reproduces much of the observed 
internal tide present in the inner basin. Several complications were 
encountered with respect to the geometrical simplifications of the model; 
however, it did substantiate that surface tides interacting with the sill at 
Sproat Narrows create internal waves of large amplitude. 


Development of meaningful criteria for dumping of potentially harmful 
substances has been approached through a series of elutriate tests, release 
experiments, and extraction experiments to determine the rates and mechanisms 
of liberation of various metals from sediments. Contaminated sediments from 
Victoria Harbour and False Creek were used. The elutriate test is based on a 
water to sediment ratio of 4:1 while the release experiments use a much lower 
concentration. The release behaviour of iron is critically controlled by 
the'dissolved oxygen, but this does not appear to be the case in the release 
pattern of the other elements. Elutriate tests on manganese, iron and zinc 
involved transfers to the dissolved phase of a very small percentage of the 
total amount of each element available in the sediment. These tests did not 
indicate, however, the amount of release which might take place at lower 
particulate concentrations in the environment. 


The selective extraction scheme operates through six regimes: interstitial 
water, exchangeable, reducible, oxidizable, moderately reducible and residual. 
The two most important reservoirsin the sediments for all observed elements 
would be the oxidizable and residual compartments, in that order. The major 
location for mercury was the residual phase where there was no evidence of 
transfer from the particulate to dissolved phase. It is apparent that at low 
particulate levels during the release experiment, metals do become mobilized 
from the oxidizable compartment. Lead, copper and mercury appear not to 
mobilize from the sediment. Although cadmium demonstrates high mobilization, 
the water concentrations are likely not to venture much above background. 


Earlier attempts have been made to use species of holothurians as a 
possible monitoring organism for heavy metals, mainly because of their size and 
abundance. Speciments of the holothurian, Molpadta intermedia, sampled from 
the Point Grey dumpsite near Vancouver did not provide adequate material to 
evaluate the utility of this species as a monitoring organism, at least at 
the analytical precision available. What evidence is available may indicate 
that there is no measurable effect of heavy metals on the benthos at the 
dumpsite. 


Sampling variability determines the number of replicates required to 
characterize a material to be dumped. In continuation of work begun last year, 
determinations have been made of the sample to sample variability in 
polychlorinated biphenyls in dredged materials of Victoria Harbour. The 
arochlors measured in samples taken during the pre-dredge survey were, in all 
cases, higher than those measured for samples removed from the barges. Grab 


samples obtained as part of the R.O.D.A.C. permit application correctly 
identified the levels of PCB present within the sampling and analytical 
error. 


A study on lead methylation has been conducted on sediments from various 
locations on the B.C. coast. The results indicate that the discharge of 
lead-bearing particulates may not result in the mobilization of lead by 
biomethylation processes under anoxic conditions. 


Future ocean dumping research will concentrate on two specific topics: 


1. Continued efforts toward the development of meaningful criteria 

for maximum permissible concentrations of contaminants in dumped materials. 
This information could be used eventually for amendment of regulations 
under the Ocean Dumping Control Act. 


2. Projects which will provide basic information to assist in the 
selection or relocation of dumpsites. 


A list of participants in the workshop is given in Appendix I. Contracts 
for 1978-1979 are summarized in Appendix II. A complete list of the contracted 
studies for the coming year, 1979-1980 is not available. 


Some of the references cited here are manuscripts produced under contract, 
and are not available for distribution. These reports as completed for the 
1978-1979 program, have been listed in Appendix IV. For further information 
on any of the Institute of Ocean Sciences series of Contract Reports 
contact the listed Scientific Authority. The locations of the Scientific 
Authorities and a key to the research institution abbreviations can be found 
in Appendix III. Pacific Marine Science Reports (PMSR) are available for 
limited distribution from the Institute of Ocean Sciences, Patricia Bay. 


This report was prepared by Dobrocky SEATECH Limited under contract to 
the Institute of Ocean Sciences Ref. DSS File No. 07SB.FP833-8-2465. 


II. OBSERVATIONS OF TEMPORAL CHANGES IN BIOLOGICAL COMMUNITIES IN ALBERNI INLET 


C.D. Levings 


Wr, Liv 
Contractor: Dobrocky SEATECH Limited 


This section covers two contracts dealing with the biological communities 
in Alberni Inlet. The major contract studying the benthic communities has 
been on-going for several years and is discussed first. The second and smaller 
contract acts as a literature search on infaunal life histories to supplement 
the data collected during the major contract. The results of this second 
contract are presented in the final paragraphs of this section. 


Summary and Conclusions 


Biological sampling in Alberni Inlet during the 1978-1979 fiscal year 
completed a project to document colonization of the dumpsite during and after 
the dumping season. This program, begun in November 1977 and conducted monthly 
until November 1978, was performed under contract by Beak Consultants Limited, 
and then, by Dobrocky SEATECH Limited (Byers and Brinkhurst, 1979). 


Two stations were monitored throughout the study. The dumpsite station 
was located on the west shore of the upper basin of the Inlet about one mile 
south of Stamp Point at approximately 56 m depth; and the "control" station 
was in about 66 m depth, about 400 m south of the dumpsite. Triplicate benthic 
grabs were taken at each of the two stations using a 0.1 m* Smith-McIntyre grab. 
Sediment from each of the six replicates was analyzed for grain-size, loss-on 
ignition, and carbon/nitrogen ratios. Measurements of water characteristics 
(temperature, salinity, and dissolved oxygen) were taken 1 m from the bottom. 
Whenever bottom conditions permitted, oxygen comsumption of the benthos was 
measured with a tethered in sttu respirometer. 


Preliminary results of temporal differences in the fauna are shown in Table 
1. Only mean values are indicated. These results are subject to revision 
until statistical analyses are complete. Conclusions regarding continued 
monitoring of the dumpsite in the Inlet, specifically, how many times per year 
and during what seasons, will be more apparent when the analyses are completed. 


Literature Review of Some Infaunal Life Histories 


Summary and Conclusions 


A contract was let to Dobrocky SEATECH Limited for the purpose of 
examining the life histories of species featured in the statistical separation 
of four "biotic provinces" in the 1976 dumpsite survey (Brinkhurst, 1979). 
There were few available data in the literature on features such as life span, 
reproductive periodicity, and fecundity but the feeding strategy and mode 
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of movement had been examined in most of the 13 species. Table 2 presents a 
summary of the characteristics. 


The "dumpsite province" was utilized by two species of polychaetes, namely 
Capttella capttata Fabricius and Sehistomeringos longicornis (Ehlers). Both 
exhibit life history characteristics enabling rapid colonization and use of 
the degraded habitat. Such characteristics include relatively high fecundity 
(>1000 ova), planktonic dispersal of larvae, deposit feeding strategy, and 
known tolerance to hydrogen sulfide. Prominent species in the "wood debris 
province" were those that used wood for food and habitat, namely the mollusc 
Xylophaga washingtona and the isopod Limnorta lignorum, both wood-burrowing 
organisms. The species from the other "provinces" were relatively unknown 
biologically and hence generalizations were not possible. 


GROUP DISCUSSION 


Rob Macdonald inquired about the error bars involved with the respirometer 
results. Depending on weather and the condition of the substrate, an 
estimation of +20-50% error was indicated: the mean deviation, though, was 
difficult to assess. The question was asked as to the comparability of the 
results from Alberni Inlet with a shallower area like Crofton (on the east 
side of Vancouver Island) and whether the two areas would demonstrate a 
similarity in species and/or statistics (Harbo, Anderson). If at all 
comparable, these 2 dumpsites might demonstrate a recolonization by similar 
opportunistic species but it was added that it was not unusual for a 
different arm in the same inlet to have a completely different community. 


Levings was asked if, after the last three years of study, he had developed 
any strong feelings concerning the impact ocean dumping was having on the 
infauna at the head of the inlet (Hoos). His response was that the dumpsite 
should remain where it is, subject to monitoring. According to several 
studies on the oxygen budget for the deeper waters of the inner basin, the 
dumped wood waste is not a large sink of dissolved oxygen compared to the 
natural reservoirs. However, this needs confirmation by refining the terms 
in the oxygen budget. 
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III. DEVELOPMENT OF MEANINGFUL CRITERIA FOR 
OCEAN DUMPING 


R.W. Macdonald 
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Contractors: Can Test Limited 
Seakem Oceanography Limited 


Introduction 


Monitoring for ocean dumping has been carried out both here and in other 
countries particularly since the "Convention on the Prevention of Marine 
Pollution by Dumping of Wastes and Other Matter" was approved in London, 
England in 1972. In that convention some substances were prohibited from 
being dumped, except as trace contaminants in other materials. The idea 
behind monitoring has been to prevent the dumping of "harmful" substances into 
the world oceans. What constitutes harmful and how to measure it has remained 
a rather intractable problem for many substances, including metals. 
Traditionally, metals in sediments have been determined by an acid digestion 
of a dried sample, this method being relatively simple, reproducible and 
inexpensive. Unfortunately this "total" metal content is very difficult to 
relate to possible environmental effects (Windom, 1975). It is, however, in 
continuing use and indeed, Canadian Ocean Dumping Regulations, where they 
specifically mention metals, refer to restrictive levels of total content 
above which a sediment may not be dumped. 


More recently this obvious problem has been addressed by the use of an 
elutriate test (Keely & Engler, 1974). Although perhaps useful for estimating 
the effects of suction dredges, it does have serious limitations, some of 
which will be brought out later. Ideally, an inexpensive analytical technique 
is required which can identify the mobile fraction of metals in sediments, and 
recently Brannon et al. (1976) have investigated the use of selective chemical 
extraction to classify sediments. Basically, this approach compartmentalizes 
metals according to the chemical treatment required to liberate them from 
sediments. It has been emphasized that it is a purely operational technique 
and therefore depends on exactly how a particular sediment responds to a 
chemical treatment. Interpretation in terms of the form of metal binding 
must remain rather vague. 


For our criteria study we have used the basic approach of Brannon et al. 
(1976) to examine the metal distribution in two local sediments, Victoria 
Harbour (VH) and False Creek (FC). We have also conducted 3-day release 
experiments in sea water with suspended sediments at approximately the 200 mg 
L7! level. After the release experiment, the sediments were retrieved and 
again subjected to the selective extraction procedure, in order to see what 
internal rearrangements the metals had undergone during the release experiment. 


Extraction Scheme 


Extraction schemes for sediment have been available for some time, 
(Gupta & Chen, 1975, Nessenbaum, 1972, Chester and Messiha-Hanna, 1970). 
It was decided to follow the form developed by Brannon etal., (1976) since 
they had made the most extensive tests and their results were directly 
applicable to the ocean dumping problem. Their scheme operates in the 
following manner: 


1. INTERSTITIAL WATER (IW) 


Care was taken to maintain anoxic conditions. 


2. EXCHANGEABLE (E) 


Metals put into solution by 1 N ammonium acetate at the surface pH 

of the sediment. Solution is defined here as passing a 0.45 y filter. 
In theory, metals located here are absorbed on organic or inorganic 
surfaces and will be mobilized in response to salinity changes. 
Release from this type of site should be rapid. 


3. EASILY REDUCIBLE (ER) 


Metals put into solution by a 0.1 M hydroxylamine hydrochloride - 
0.01 M hydrochloric acid solution. This phase consists of oxides, 
hydroxides and hydrous oxides of manganese and iron in the form of 
coatings or micro particles. Scavenged trace elements and phosphate 
may also be present (Brannon etal., 1976). 


4. OXIDIZABLE (0) 


Metals put into solution by a 95°C digestion with 30% hydrogen 
peroxide at a pH of 2.5. These metals originate largely from 
. organic complexes and sulphides. 


5. MODERATELY REDUCIBLE (MR) 


Metals put into solution by a sodium citrate - sodium dithionite 
solution. Metals released during this step are probably associated 
with hydrouns iron oxides which were not reduced during the "easily 
reducible" chemical treatment. 


6. RESIDUAL (R) 


Essentially the metal that remains after the above extractions and 
determined by digestion with hot aqua regia. Metals located here 
tend to be bound in the matrix and so generally are not released 


except under extremely harsh conditions not normally encountered in 
the environment. 


In addition to this selective extraction procedure, the sediments were 
subjected to several other determinations, including % organic carbon, size 
distribution, free and total sulphide, cation exchange capacity, pH and Eh. 
A total analysis for metals was also performed on’ a subsample of each 


homogenate in order to estimate the mass balance of the selective extraction 
scheme. 


Results 


Comparison between elutriate tests, release experiments and selective 
extraction is not possible unless some common base iis used. For all of the 
metals, the data have been put into concentration units of ug g-! (dry weight). 
In that way, releases expressed in these units can be viewed directly in terms 
of the amounts of metal being transferred to the liquid phase during the 
elutriate test and alsoto the amounts of metals residing in each of the 
compartments as indicated by the selective extraction procedure. In actual 
fact, different solid to liquid ratios have been used: the elutriate test is 
based on a 4:1 (V:V) water to sediment ratio, whereas the release experiments 
performed here have used 40g (wet weight) of sediment in 200L of sea water. 
Nevertheless fluxes and concentrations expressed in terms of ug g! dry 
weight sediment sediment should be directly comparable. 


Lron 


The results for iron are summarized in Figures 1 and 2. In referring to these 
it will be noted that for both Victoria Harbour (VH) and False Creek (FC) 
sediments, run #1 and run #2, which were supposed to be duplicate experiments, 
do not compare well. The reason for this is both revealing and simple. The 
cores, when collected, were kept sealed, and when delivered to Can Test were 
maintained under Ny in a glove box, a procedure known to be essential for 
meaningful pore water measurements (Poldski & Glass, 1974). At this time 
subsamples were removed for various chemical tests and also a subsample was 
placed into a sealed polyethylene container and sent to Seakem for the release 
experiments. Anoxic integrity was maintained in this manner. At Seakem, 
subsamples were removed for run #1, weighed and added to the stirred 200L 
container. The remaining material was resealed in the polyethylene container 
and stored at 4°C in the dark for a period of almost 2 weeks. Since the 
bottles were opened under a normal atmosphere, anoxic integrity was lost at 
this point. At the start of run #2, therefore, the sediments had been 

exposed to oxygen for a two week period. 


In Figures ] and 2 it can be seen that, for run #1, there is an initial 
period of 1-2 hours during which very little release takes place, and perhaps 
there is even a small amount of scavenging. With 4 hours, noticeable 
releases of iron took place for both sediments, with a maximum release 
ocurring after about one day. But within three days, the released iron had 
been scavenged and the dissolved levels returned to near background. Sediments 
which had already been exposed to oxic conditions for two weeks (run #2) 
behaved differently, initially releasing only a small amount of iron. The 
reproducibility of this effect between FC and VH sediments is excellent. 
Clearly, much of the process ocurring in the water column during run #1 was 
also ocurring in the material stored in the bottle and destined for run #2. 


Of particular importance are the elutriate test results. These give no 
information on the longer term iron releases from anoxic sediments, and one 
might have predicted this since the test involves shaking the water-sediment 
mixture for only 4% hour. Furthermore the transfer of iron from the solid to 
the dissolved phase (ug Sooke dry weight units) is very dependent on the 
concentration of particulate material. Duriug the elutriate test only 0.81 
ug, g-? and 3.8 ug gl was put into the water for the VH and FC sediment, 
respectively. Looking at Figures 1 and 2 it can be seen that these levels of 
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transfer hardly dip at all into the iron reservoir in the sediment. On the 
other hand, during the release experiment (run #1) as much as 1000 ug go) 
transferred out of the solid phase, or approximately 4% of the iron reservoir 
was being affected. 


Examination of the selective extraction data shows the comparison of 
before and after partitioning of iron. Releases occurring during run #1 
cannot be derived from interstitial water (IW) or from the exchangeable portion 
(E). In both sediments, the easily reducible (ER) portion comes closer to the 
amount in question, and indeed for both sediments there is a significant 
depletion of iron from this compartment after the experiment. Analytical errors 
probably do not allow one to make any definitive statements about the remaining 
compartments since inter-core duplication was not observed. It should be 
noted that the bulk of the iron in both sediments is found in the oxidizable 
phase (sulphides and organics) and in the residence material. 


Manganese 


The behaviour of manganese is quite different from that of iron (Figures 
3 and 4). In both the VH and FC sediments there was an immediate (1 hour) 
loss of dissolved manganese from the water. This was followed by a release 
from the sediments back into the dissolved state which continued throughout 
the rest of the experiment. Run #2 (keeping in mind the exposure of the 
sediments to oxygen) did not show the initial scavenging from the water, but 
the trend thereafter was for dissolved manganese to increase. In run #2 of 
the VH release experiments, the last data point (3 days) shows that the 
dissolved manganese had drastically decreased and was, in fact, below the 
initial background level. Whether or not this is genuine or an experimental 
error is difficult to say since it was not repeated in the other release 
experiments. What is interesting is that during three of the release 
experiments, the dissolved manganese levels approached those occurring initially 
in the interstitial water (21-23 mg L~!). For comparison, the final 
concentrations of the four release experiments were 19.2, 7.6, 19.9 and 21.1 
ng L-!, To further emphasize the importance of this sediment: water ratio 
the elutriate test gave final dissolved manganese concentrations of 53 and 39 
mg L7!, Actual losses from solid to liquid during the elutriate test were, of 
course, remarkably small being only 0.46 ug g-! and 0.20 ug go! . The elutriate 
test therefore may have approximately reflected the water concentrations of 
manganese expected, but a detailed knowledge of the amount of manganese 
released cannot be known without information on the time - concentration path 
of the suspended sediments. 


The selective extraction experiments showed that the bulk of the manganese 
was in the oxidizable portion and the residual portion of the sediments. 
Manganese concentrations in the interstitial water and in the elutriate test 
were generally low, and referring to Figures 3 and 4, the transfer of manganese 
to the water during the elutriate test must also be considered low. The 
exchangeable portion could have contributed sufficient manganese to produce 
the releases observed during the elutriate test. The release experiments 
cannot be explained so easily since at the end of the experiment a significant 
fraction of the manganese was in the dissolved phase. It would appear that 
portions of the oxidizable and/or residual manganese must have been mobilized 
during the three day release experiment. 


Zine 


The behaviour of zinc appears to be more closely allied with that of 
manganese than that of iron. The similar trend of slow increase throughout 
the release experiment is evident in Figures 5 and 6. As was the case with 
manganese, the zinc concentration tended to approach a common level of 6-8 
Lg L7! in all four experiments. The results of the elutriate tests showed 
that more zinc was mobilized from the VH sediments than from the FC ones, 
and the reason for this may be that the VH sediments released the zinc faster. 
It is clear, however, that not much zinc is being transferred to the dissolved 
phase. ; 


The majority of the zinc in both the VH and FC sediment appeared in the 
oxidizable portion with somewhat lesser amounts in the residual phase. Zinc 
mobilization during the release experiment can only be explained by 
proposing that much of the oxidizable portion of the zinc reservoir was being 
mobilized. Both mass balance during the release experiment, and comparison of 
the "before" and "after" distribution of zinc in both sediments bear this out. 
There is even evidence that some zinc is transferred from the oxidizable phase 
into the exchangeable and easily reducible phases of the sediment during the 
three day residence in oxic waters. For the elutriate tests the amounts of 
zinc being mobilized are so low that one does not require loss of the zinc 
from the oxidizable phase to explain the results. However it is 
interesting tonote that some correlation (r=0.516) between the oxidizable phase 
and elutriate test results was observed by Brannon et al. (1976) based on 
several cores from different locations. The present results seem to give some 
very direct evidence that leaching of zinc from the oxidizable portion can take 
place. 


Copper 


The results of the selective extraction measurements for this element. are 
presented in Table 3. The total levels of copper in VH sediments appears 
quite high, 196 ug g-!, but falls into the range determined during the study 
(Byers and Brinkhurst, 1979). The copper levels are about half this in the 
FC core but in both cores the majority of material is located in the oxidizable 
phase with most of the remaining copper being accounted for by the residual 
phase. 


The release experiments did not show any measureable amounts of copper 
being mobilized and,in fact, small amounts may even have been scavenged. 
Dissolved copper levels remained at or near detection levels after the sediment 
was mixed into the sea water. During the elutriate test detectable but very 
small amounts of copper appeared to be released with the scale of release 
being about 0.02 - 0.04 ug g 


Cadmium 
Table 3 includes the selective extraction data for cadmium and it is 


apparent that, like copper, the majority of the metal is found in the oxidizable 
portion. Due to the low levels (below detection limits) it was impossible to 


assign the remaining cadmium more accurately, but the determination of total 
cadmium gave levels of about 1-1.2 ug gt. enough to prohibit dumping of the 


material. 


The data obtained for cadmium during the release experiment were not 
good enough to make a case for either release from or adsorption onto 
suspended sediment. Detection limits and contamination are two problems which 
make the task of measuring dissolved cadmium difficult, if not impossible. 
However, it should be noted that the amount of cadmium remaining in the 
oxidizable phase after the release experiments was somewhat less than that 
at the beginning, a behaviour similar to that observed Tou ofits. UMEGOLS 
basis, the evidence favours some loss of cadmium from the solid phase. The 
elutriate test gave no evidence of cadmium going into solution, even from these 
sediments which were relatively enriched in cadmium. In fact, all cadmium was 
below detection limits. 


Lead 


Selective extraction data for this element is also found in Table 3 where 
it can be seen that the lead levels in VH sediments are high but fall into 
line with previous determinations (Byers and Brinkhurst, 1979). Following the 
behaviour of copper, zinc and cadmium, the majority of lead was found in the 
oxidizable portion of the sediment. Most of the remaining lead was accounted 
for by the residual phase. 


It is very difficult to demonstrate either release or adsorption of lead 
during the release experiment as was true of Cadmium, and lead remained below 
the ug L7-! level at all times. In any event it would seem that no more than 
1% of the total lead was mobilized. 


Arsenic 


Results from the selective extraction procedure are presented in Table 3 
and again it is apparent that levels are often at or below detection limits, 
making mass balances difficult to perform. Most of the arsenic appears in the 
oxidizable phase, and the residual compartment is of secondary importance. 

The release experiments did not yield any useful data for this element, 

but on the basis of the "before" and "after" selective extraction procedure, 
almost half of this element tied up in the oxidizable phase was lost, while the 
residual phase remained untouched. The behaviour of arsenic therefore seems 

to correlate with metals described above. 


Mercury 


Table 3 summarizes the data for this element which behaved somewhat 
differently from the other metals. The majority of this element was tied up in 
the residual phase and in keeping with this there was no evidence of transfer 
from the particulate to dissolved phase during the release experiment. In 
fact dissolved Hg levels tended always to drop slightly after the sediment 
was added to the sea water. In addition the elutriate test did not reveal any 
mercury mobilization. 
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Conclusions 


Not all of the data obtained have been presented since this report is 
intended to present concisely the salient features, particularly as they apply 
to present methods of ocean dumping monitoring. 


The first important point is the effect of storing sediments in contact 
with oxygen. Obviously, as Figures 1 and 2.clearly show, the behaviour of 
iron is controlled by this. However no clear evidence was found for any 
change in the release pattern of the other elements investigated,at least 
within the detection limits and precision involved in the respective 
measurements. 


In view of the behaviour of manganese, iron and zinc, the elutriate test 
did not seem to be indicative of the amount of release which might take place 
at lower particulate concentrations in the environment. Transfers to the 
dissolved phase during the elutriate test generally involved very small 
percentages of the total amount of each element tied up in the sediment. One 
might therefore erroneously assume that much of the element in question was 
unavailable. An accurate picture of the time-concentration history of dredge 
spoils after dumping would be required before one could calculate the release 
budget. Maximum concentrations however might be quite successfully estimated 
by using the elutriate test, although the half-hour time period may prove to 
be rather arbitrary. 


For all of the observed elements the two most important reservoirs in 
sediments would be the oxidizable and the residual compartments, generally 
in that order, except for manganese and perhaps iron. Mercury was different 
in that the major location for this element was the residual phase. 
Differences in concentration before and after the release experiment indicate 
that, at low particulate levels during the release experiment, metals do 
become mobilized from the oxidizable compartment. Actual amounts of metals 
found in the dissolved phase during the release experiment also indicate that 
the exchangeable and interstitial water phases do not contain enough of the 
elements to explain the results. In the case of lead and copper, mobilization 
does not appear to occur, or if it does it includes only a small percentage of 
the total metal contained in the sediment. Mercury also does not appear to 
be mobilized, in agreement with the fact that the majority of this element was 
found in the residual phase. Cadmium does, however, appear to be mobilized to 
an extent that perhaps 50% was lost during the three day experiment. Both 
the elutriate and release experiment though, indicated that water concentrations 
of cadmium are not likely to venture much above background. 


Group Discussion 


The results of the experiments identified Cd mobilization in sediments 
as being of possible concern. With respect to this, two attendees (Thompson, 
Waldichuk) asked if the elutriate tests conducted could be extrapolated to 
determine uptake from the interstitial water by animals. Apparently the tests 
could not, since such conditions would involve a totally different environment 
more closely linked with sedimentological work. Harbo asked how sediments 
containing a large amount of wood debris would compare to those sediments used 
for the experiments. The False Creek sediments, in fact, did contain wood 
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debris and the general tendency was towards mobilization of the metals. Since 
there could be a short-term depression of 02 in the water colum during and 
subsequent to dumping, Levings asked about the time frame of biological 
availability (concommitant with the release pattern of Fe in oxic conditions). 
After a typical dump, the metals would be available to planktonic organisms 
while the material was in the water column and could be released from these 
sediments up to 3 days afterwards; however, once the dumped material struck 
bottom a different biological and chemical process would be involved in 
reaching an equilibrium. What test, then, could be recommended as a means of 
determining biological uptake of these metals (Waldichuk)? A weak acid leach 
was suggested as a better method than the standard elutriate test for making 
the metals available (although the Eh and pH must be closely monitored). Young 
asked about the influence of reduced salinity such as that found in estuaries 
(eg. 5%) » Macdonald acknowledged that a change in salinity could cause an 
alteration in release times and a possibility of rescavenging of the element 
back to the solid phase. For example, Zn and Cu have been found to mobilize 
when the fresh water of the Fraser River encounters the salt wedge. 
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IV. HEAVY METALS IN Molpadia tntermedia FROM 
THE POINT GREY DUMPSITE 


J.A.J. Thompson 
LL. 0.8. 


Contractor: Can Test Limited 


Introduction 


In previous reports (Brinkhurst, 1978; Thompson and Paton, 1978) the use 
of the holothurian, Molpadia intermedia as a possible monitoring organism 
for heavy metals in the Point Grey Dumpsite area has been described. The 
previous study produced inconclusive results. Reasons for this were basically 
due to the lack of a uniform statistical sampling of the dumpsite, lack of 
sufficient samples, and, the questionable precision of the analytical data 
provided by the contractor. 


In July 1977 a cruise to the dumpsite on the C.S.S. Vector provided 
samples of M. tntermedia from the 17 stations shown in Figure 7. At least 5 
specimens were collected from each station. Controls were obtained at station 
F off the Sechelt Peninsula (49921.6' N, 123°34.8' W). 


Sampler casts at McCall Bank (49920.4' N, 123935' W) and Satellite Channel 
(48942.5' N, 123930' W) were unrewarding. Control comparability was seriously 
hampered for this reason. 


All samples were obtained using a Smith-McIntyre grab (Kahlsico, San 
Diego). Sediments were passed through a 1.0 mm nylon mesh screen. Specimens 
were vigorously scrubbed to remove adhering sediment, weighed, and immediately 
after they were dissected in a clean-air hood. The longitudinal muscles 
were removed from each specimen and frozen in separate acid-washed vials. 
Preparatory to analysis, all samples were freeze-dried. The contractor was 
provided with a total of 83 coded samples, including NBS standard bovine liver 
and orchard leaves as analytical controls. Samples were digested in aqua 
regia and analyzed by flame or flameless atomic adsorption methods as required 
for Cu, Zn and Cd. 


Results 


Analytical results are summarized in the following table: 


Range * Control 
(ug g ! dry weight) (ug g7! dry weight) 
Copper & le ae 9.1 6.6 
Zine 121.0 - 146.0 129.0 
Cadium 0.4 - 4.4 1.9 


a 
Range indicates 5 or more determinations per station. 
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As was the case in the first study little can be concluded about trends 
in concentration in the dumpsite area for any of the metals studied. It 
should be noted that Cr and Pb were not determined in this present study 
because of financial restraints. These elements also provided the least 
precise data in the previous work. 


Copper does exhibit a slight trend to increasing concentrations in the 
western sector; however, it is not significant. Moreover, the control lies 
midway in the range of dumpsite specimens. 


An interesting aspect of the copper data is that the means are considerably 
lower than those obtained in the 1977 study. Overall copper mean for this 
work is 5,99 2.1775 vg go! dry weight against 26.0 + 14.0 for the previous 
data. This significant difference may reflect greater care taken in sample 
preparation and/or greater analytical accuracy and precision. There is no 
a prtort ecological reason for this decrease. 


Zinc also exhibits a similar significant drop in the overall mean from 
ile te oo 0 tO) 37s 0Cfe 0 UE gi dry weight for 1977 and 1978, respectively. 
Here again the increase in precision may reflect the greater care in sample 
preparation and analysis. 


Zinc is the only one of the three elements that exhibits noticeably 
higher data for the dumpsite when compared to the unfortunately scant control 
data. It is, however, not possible to indulge in a valid statistical exercise. 


Data for cadmium are consistent with those obtained previously. Overall 
means for this and the earlier work are 1.4 + 1.0 and 1.7 + 1.4 ug on diy 
weight, respectively. 


Conclusions 


Although analytical precision has improved it is still not possible to 
observe differences (if any) outside of the "noise" or background level of 
natural variability and analytical error. We were limited by the size of 
samples which were smaller than would be desired. 


The holothurian M. intermedia does not fulfill expectations as a possible 
monitoring organism for several reasons; 

1. it is insensitive to changes in sedimentary metal concentrations; or 

2. there are no significant increases in sedimentary levels; or 


3. it does not provide enough tissues from which valid data can be 
obtained. 


If the first or third conclusions are true then it might be concluded 
further that another organism should be tested. An alternative observed in 
sufficient numbers was the sea urchin Brisaster sp. 


If the second conclusion were correct it would be assumed that there is 
no measurable effect upon the benthos from dumping activities. 
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Group Discussion 


Waldichuk mentioned that some of the cadmium levels in Molpadia intermedia 
were above those allowed by the Food and Drug Directorate for human consumption. 
This particular sea cucumber, however, was not expected to reach market value! 
In response to a question from Levings, Thompson mentioned that levels of 
copper had decreased between this study and the previous study done by E.P.S. 
where results from their routine sediment survey at Point Grey in 1978 
(analyzed by Fletcher's techniques at the University of British Columbia) 
also indicated that metal values have all shown a reduction. 
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V. VARIABILITY OF CONTAMINANT CONTENT IN DREDGED SPOILS 


R.W. Macdonald 


r.0.8 


Contributors: C.S. Wong and P. Gregory 
Contractor: Seakem Oceanography Limited 


Introduction 


During the last fiscal year, an examination of the reliability of sampling 
for ocean dumping monitoring was performed. As has already been reported 
(Byers and Brinkhurst, 1979) the sampling strategy involved obtaining seven 
grab samples at random locations throughout the dredge site before the 
commencement of dredging. In addition, four of the barges were sampled at 
seven randomly selected locations within each barge after it had been filled 
with spoil. Each sample was then split into three subsamples, stored in 
containers appropriate to the analysis to be performed and then randomly 
numbered so that the analyst would be unaware of replicates. This year the 
project was completed by examining the levels of polychlorinated biphenyls 
(PCB) in the 105 samples. 


PCBs are a family of compounds consisting of over 200 individuals each 
being a particular chlorinated biphenyl generated by replacing the hydrogen 
at any combination of the positions numbered 2-6 and 2'-6' by chlorine. 
PCBs produced industrially are generally a mixture of many different PCB 
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compounds and have been designated by referring to them as arochlors. A 
four figure number is added to indicate the number of carbon atoms involved 
and the percentage replacement of hydrogens by chlorine. For example, archlor 
1254 contains 12 carbons (' biphenyl) and has an average 54% replacement of the 
hydrogens. Quantitation of biphenyls is therefore difficult because, in 
principle, a single number is produced (ug go!) which covers a host of different 
compounds. Quantitation has been performed by using a set standard, arochlor 
1242 for example, which comes closest to the distribution and range of 
compounds in the sample at hand (Harding et al-,1978). Although quantitation 
is not strictly accurate (it is like giving a grand number for the concentration 
of 100 different compounds each of which has not been calibrated), good 
precision can be achieved with a decent match between standard and sample. 
For the present results, quantitation was carried out by dividing the 
chromatogram into three areas each of which closely matched a particular 
arochlor. In addition decachlorobiphenyl was also quantitated and all of 
the results are summarized in Table 4. Levels of PCB determined last year by 
Willis, Cunliffe, Tait and Company Limited for E.P.S. permits (4443-0375) as 
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part of the normal R.O.D.A.C. permit process, are also included and show rather 
respectable agreement with the totals summarized here. 


The 105 samples were decoded and subjected to analysis of variance, the 
results of which are presented in Table 5. For all of the arochlors, and the 
total, the between-replicate variability accounts for over 40% of the variation 
observed. The variation within each barge accounts for approximately a further 
25% with the remaining variation being due to differences between barges. The 
ANOVA results show that variability within barges is significant only for the 
DCB fraction but differences between barges are real for all arochlors at the 
99.9% confidence level. A point worth noting is that the arochlors measured 
in samples taken during the pre-dredge survey were, in all cases, higher than 
those measured for samples removed from the barges. 


Conclusions 


Due to the large number of different PCB compounds, exact quantification 
would be prohibitively expensive, and probably not necessary for ocean dumping 
monitoring unless a particular PCB compound was found to cause exceptional 
environmental problems. The ANOVA results show that there are genuine inter- 
barge differences but an examination of the data showed that in all cases the 
pre-dredge sampling programme provided samples with higher PCB levels than 
those found in the barges. Grab samples obtained as part of the R.O.D.A.C. 
permit application correctly identified the levels of PCB present within the 
sampling and analytical error. 


Group Discussion 


In response to a question from Levings, Hoos pointed out that PCB levels 
“in Victoria Harbour were much lower than those reported from Puget Sound 
or Porpoise Harbour. Macdonald concluded the discussion by emphasizing 
the fact that the study on Victoria Harbour and False Creek sediments (PCB 
analyses completed this year and the trace elements and hydrocarbon analyses 
completed last year, 1978) demonstrated that adequate sampling done prior to 
dredging in a new area could reduce the monitoring required subsequently at 
the dumpsite. 
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VI. METHYLATION OF LEAD IN MARINE SEDIMENTS 
J.A.J. Thompson 
Contractor: Beak Consultants Limited 


Introduction 


It is known that a number of other elements besides mercury can form stable 
methylated compounds through biological activity in natural media. The 
mechanisms for the production of CH3Hgt and (CH3)) Hg are understood to some 
extent; however, little is known about the mechanisms and importance for the 
methylation of selenium, lead, arsenic and tin (among others). 


Studies by Chau and coworkers at the National Water Research Institute in 
Burlington (Wong et al.,1975; Chau and Wong,1978), Jarvie in England and Huber 
in Germany have concentrated on the methylation of lead in lacustrine, and 
fluviatile sediments as well as in axenic bacterial cultures derived from them. 


Contradicting results have been reported by Wong et al.(1975) and Jarvie 
et al.(1975). Recent work has clarified the situation to some extent. It is 
now known that both Pb(II) and Pb(IV) species can be methylated bacterially. 
Lead(II) methylation proceeds with difficulty, probably because an oxidative 
step is required. R3Pb cations,on the other hand,undergo facile methylation. 
Jarvie et al.(1975) first suggested that this reaction was abiotic (i.e. purely 
chemical). Studies by Chau et al. (1978) and Huber et al. (1978), however, 
clearly show that both biotic and abiotic reactions are extant with the latter 
being the dominant process. 


Little is known about the importance of methylation of lead in marine 
ecosystems. Because total lead concentrations in sea water are extremely 
low (5-15 ng L-!) it would not be expected (nor could it be easily determined) 
that methylation plays an important role in the water column. Lead in natural 
sediments is in a predominantly immobile phase and here again it might be 
predicted that little activity would be expected. (See pp. 12 of this report 
to identify mobility of lead). 


Our interest in the lead methylation capability of sediments arose out 
of the studies on metals release from submerged mine tailings deposits in 
coastal inlets. We had found previously that ambient Hg in sediments is 
available for methylation, and were interested in whether lead may exhibit 
similar properties. 


For this study, sediment samples from four locations on the B.C. coast 
were collected: 

1. Burrard Inlet (English Bay). 

2. Goletas Channel (Northern Vancouver Island). 


3. Alice Arm (site of inactive molybdenum mine with high lead in 
tailings). 
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4. Granby Bay (site of abandoned copper mine at Anyox). 
Methods 


Collection was made with a Smith-McIntyre grab. A subsample from each 
grab was placed into a dark green plastic bag, purged of air, and sealed 
tightly. They were kept at 4°C until needed for the experimental program. 


A series of flasks containing sediment from each of the four areas was 
incubated anaerobically at 15° for periods up to 4 weeks. Sampling from the 
head-space gases was made with a septum at 0, 1, 7, 72, 168, 336 and 672 hours 
into the experiment. Three groups of 4 flasks were spiked with trimethyllead 
acetate, lead(II) nitrate and lead(II) acetate, respectively. A fourth set 
of 4 flasks served as controls. Each of the four flasks in all four groups 
contained solutions according to the following scheme: 


1, 75 ml of unsterilized sea water 

2. 75 ml of sterilized sea water 

3. 75 ml of sea water, 1.0 ml of nutrient broth and 1.0 ml glucose 
4. 75 ml of (3) with sterile sea water substituted. 


Analysis for (CH3), Pb in head space gases was performed using a GLC 
with an electron capture detector. Confirmation of the presence of (CH3), Pb 
(TML) was made using a Finnigan GC/MS. Detection limit for the GLC method 
was 0.23 ug 1 ie 


All samples were also analyzed for total lead (aqua regia and nitric 
acid extractable), total organic carbon, moisture and total sulphide. 


Results and Conclusions 


Results of the incubation experiments can be summarized as follows: 


1. No TML was detected in the controls. This suggests that lead in 
sediments is unavailable to short-term methylation processes. 


2. TML was produced readily in the presence of (CH3) 3Pbr. Over 90% 
of the available Pb was present in the head space as TML at the end 
of the experiment under all four conditions. 


3. In the presence of Pb(NO3)2, TML was produced to a maximum of about 
IAG AW: L7! or 0.026% of the available lead. 


4. In the presence of Pb(OAc)9, no TML was produced above the detection 
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5. The presence of added nutrients had little effect upon methylation 
rate. Some slight differences in rates among the sediment types 
was noted. The sediment from Goletas Channel was the most productive 
in methylation of Pb(NO3); however, differences were not significant. 


6. No methylation of Pb(NO3). or (CH3)3Pbt occurred in sterilized controls. 
Addition of NajS resulted in very rapid production of TML from (CH3) 3 
Pbt but none from Pb(NO3)>. 


Results of this study indicate that if methylation of lead occurs at 
all under natural conditions it must result in the production of very little 
TL. When it is considered that incubation temperatures were about three 
times greater than those encountered in coastal sediments, any methylation 
process would proceed very slowly. 


The presence of lead in tailings-derived sediments from Alice Arm at 
levels up to 14 times those in the three other samples did not result in 
enhanced TML production. This would suggest that the lead (probably present 
as PbS) is not available to bacterial action in anoxic sediments. A recent 
report by Harrison and Laxen (1978) suggesting that tetraalkyllead in rural 
air was in part due to outgassing from intertidal mud flats, may indicate 
that oxygen exposure is required. 


With respect to ocean dumping, our results would indicate that the 
discharge of lead-bearing particulates may not result in the mobilization 
of lead by biomethylation processes under anoxic conditions. 


Group Discussion 


Considerable discussion followed on lead and its various compounds. With 
little mobilization occurring with the Pb found in mine tailings, the question 
was raised whether other forms of Pb reacted in a similar manner (Hoos). 
Thompson replied that Pb in a more available form could behave quite 
differently. For example, in these experiments, all the sediments tested 
required the addition of lead before methylation was induced. It was 
postulated that lead compounds with a greater degree of covalency in the lead- 
anion bond might permit less or no methylation. In response to a question 
from Levings, Thompson indicated that Pb used for dating cores was actually a 
Pb isotope, not stable lead. 


Waldichuk asked if this work was conclusive, but Thompson felt that 
better results could be obtained by means of a continuous flow apparatus where 
the carrier gas would collect everything that became mobilized, and in 
greater quantities. He also felt that more information was needed on the 
microbiology of methylation. 


Young asked if the bugs were the means by which methylation occurs. In 
fact, this is true; however, organisms can also work as a means of 
demethylation depending on the strain of bacterium. 


Waldichuk mentioned that the several compounds of lead have very different 
stabilities: (CH3), Pb (TML)is stable in sea water while (CH3)2 Pb is unstable. 
Also, lead oxides introduced to the sea via car exhausts are known to be 
quite non-reactive. Thompson acknowledged this and stated that Chau found 
this particular form of Pb didn't react because it was insolwle. 
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VII. AN EXAMINATION OF THE VARIABILITY OF 
UPWELLING AND ITS RELATIONSHIP TO 
THE FLUSHING OF ALBERNI INLET 


J.A. Helbig 
Contractor: Beak Consultants Limited 


General Introduction: . W.H. Bell, 1.0.5. 


The dumping of dredge spoil in the inner basin of Alberni Inlet imposes 
an oxygen demand on the deep water of the basin. The evidence available 
to date suggests that flushing or renewal of this water occurs on an annual 
basis so that there has been no long-term accumulation of deleterious effects 
due to the oxygen demand. However, the question arises as to the invariability 
of the flushing and the consequences of a lack of renewal in any particular 
year. Failure of the renewal mechanism might require suspension of dumping 
permits in the Inlet for a certain period of time to avoid additional oxygen 
depletion problems. 


The flushing of the Alberni Inlet system appears to result directly 
from upwelling on the west coast of Vancouver Island and is, therefore, 
related to the annual cycle of winds in the area. Beak Consultants Limited 
were awarded a contract to examine the problem and attempt an estimate of the 
likelihood of failure of the renewal mechanism. (Beak Report, 1979). The 
study was carried out by Dr. J. Helbig, who has summarized his results below. 
Dr. Helbig examined the relationship between coastal winds, surface sea water 
density and known renewal events in Alberni Inlet, insofar as the available 
data permitted. As a result of this examination, he suggests that upwelling 
is not so much due to the forcing of ocean circulation by north-west winds 
as it is a response to the relaxation of south-east winds. On the basis of a 
particularly strong upwelling event, which was responsible for basin water 
renewal in 1966, he establishes that there is some likelihood that an event of 
similar strength will not occur every year. However, he concludes that lesser 
events are adequate to cause renewal, in many instances, and one of sufficient 
strength for this purpose is almost certain to occur every year. 


Introduction 


That coastal upwelling might lead to the flushing of Alberni Inlet was 
first suggested by Tully (1949). Bell (1976) has clearly described how 
upwelling is related to deep water renewal and has argued on the basis of 
available, although scanty, data that the inner basin of the Inlet is replaced 
annually with oxygen "rich" water. Since local industry places a large 
biological oxygen demand on these waters, the importance of a failure in the 
presumed annual cycle is clear. This report summarizes an attempt to estimate 
the probability of a lapse in this cycle. Note that these ideas are not 
restricted to Alberni Inlet, but also apply to other inlets on the west coast 
of Vancouver Island. 


Pt 


The basic concept is deceptively simple. As the prevailing winds shift 
to southeastward in spring, the joint effect of the wind stress and Coriolis 
force is to drive surface waters away from the coast. To replace this 
water, deep (and more dense) water is upwelled from below. If the density 
of the upwelled water at sill depth exceeds that of the inlet, it will flow 
into the inlet and displace the resident water upwards. 


There are three reasons why in any given year, the deep waters of Alberni 
Inlet may not be re-oxygenated: 1) the winds might be of improper direction 
or insufficient intensity or duration to induce sufficient upwelling; 2) the 
upwelled water could be too low in density; or,3) the upwelled water may 
be deficient in dissolved oxygen. Of these three, only the first is directly 
related to upwelling, and it is the only one examined here: the others are 
concerned with the "climate" of the offshore waters. 


The available data consist of a relatively sparse set of oceanographic 
measurements taken in the two basins of Alberni Inlet and in the offshore 
zone; winds recorded at several locations on the west coast of Vancouver 
Island by the Atmospheric Environment Service; and a 41 year time series of 
surface temperature and salinity taken at the Amphitrite Point lighthouse 
station at the mouth of the Inlet. Included in the oceanographic data are the 
results of a series of 19 monthly cruises made in the inner basin in 1965-67. 


Results and Conclusions 


The wind stress calculated from winds at Tofino, the surface density 
computer from temperature and salinity fields at the Amphitrite Point 
lighthouse station, as well as the runoff from the Somass River have been 
examined to determine the relationship between coastal upwelling and renewal 
of the deep waters of Alberni Inlet. The major obstacle to the achievement 
of this goal is the severe lack of synoptic oceanographic data from Alberni, 
particularly from the outer basin which is directly connected to the Pacific 
Ocean. This deficiency of data prohibits the determination of the time and 
duration of renewal. The existing data are, however, consistent with an 
annual renewal of the deep waters of the inner basin beginning in late spring 
or summer. Obviously the outer basin is also flushed at least annually. Line 
P data clearly indicate the upwelling and downwelling of isopycnals. 


Winds tend to blow down the coast (from the northwest) during the months 
of May through August and up-coast otherwise. Local topographic features 
can alter winds measured on land, particularly near inlets. For example, 
predominately southeastward winds tend to possess a positive onshore component. 
When the winds are strong, the alongshore component of the wind stress always 
greatly exceeds the onshore component. The comparison of winds measured at 
Tofino and at Estevan Point indicates that the wind field is comparatively 
homogeneous. 


The linear correlation between the wind stress at Tofino and the surface 
density at Amphitrite Point is always small. This means that the long time 
series of surface temperature and salinity is of limited value in determining 
the relationship between upwelling and flushing. Several factors may 
contribute to the low correlation: fresh water runoff and its concomitant 
variation with the winds, local topographic modification of the winds, the 
depth of the mixed layer, non-linear effects, and wave-like disturbances that 
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are generated outside of and propagate into the region of interest. 


Visual examination of the time series of wind stress, surface density, 
and river outflow reveals distinct upwelling - downwelling events. A 
comparison of these variables with density and dissolved oxygen content data 
from the inner basin of Alberni Inlet convincingly demonstrates that upwelling 
leads to flushing. The time scale for the onset of renewal following a shift 
to upwelling-favourable winds may be a week or less. Large scale flushing of 
the inner basin appears to take at least a month. 


Based on the above comparison, estimates are made of the strength and 
duration of southeastward winds necessary to induce renewal. The probability 
for any one year period, that a 30-day interval will occur during which the 
average southeastward component of the wind stress will be of equal magnitude 
to the southeastward wind stress observed during May 1966 (when the renewal 
occurred) is about 75 percent. Since the most negative monthly and weekly 
averaged southeastward wind stress occurred during May, the calculated 
probability will tend to be underestimated and hence the likelihood of a 
lapse in the supposed annual renewal cycle is small. Based on 1965 values, 
this probability is essentially one. 


Any future observation program relating to the flushing of Alberni Inlet 
or the other inlets on the west coast of Vancouver Island should include the 
mooring of current meters near all sills and in all connecting waterways to 
the ocean. Moreover these meters should be equipped with self recording 
conductivity and temperature sensors. This program should be supplemented 
by regular synoptic cruises. In addition, provision should be made for 
the adequate measurement of the winds. 


Group Discussion 


In response to an inquiry by Hoos, regarding the regularity of renewal 
of the inner basin of Alberni Inlet, Helbig indicated that there was a 
very strong probability that an upwelling event of sufficient strength to 
initiate flushing, would occur at lease once per year. Anderson pointed out, 
however, that the probability of flushing is not independent from year to year. 
For example, the filling of the inner basin with very dense water in one year 
might reduce the probability of flushing in subsequent years. To a query from 
Levings, Helbig replied that better temporal data on density would improve 
the capability for predicting the likelihood of renewal. 


Buckley suggested that the data might be suseptible to cross spectral 
analysis as an alternative to correlation analysis. Helbig agreed, but advised 
that the scope of the project was not such as to allow it. 


Young asked if every inlet along the west coast of Vancouver Island would 
flush in synchronism. Helbig stated that this would probably happen, with the 
possible exception of Quatsino Sound. Waldichuk asked whether a large wave 
(tsunami or seismic sea wave) would influence flushing, but the effect would 
apparently be dependent on the resonance period. 
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VIII. AN EXAMINATION OF EVIDENCE FOR AN 
INTERNAL TIDE IN ALBERNI INLET 


J.R. Buckley 


Contractor: Seakem Oceanography Limited 


General Introduction: W.H. Bell, 1.0.5. 


Time series measurements of water properties made in Alberni Inlet in 
1977 revealed the existence of a large internal tide in the inner basin, a 
region which includes the designated dumpsite. The effect of such tides is 
to cause fluctuations in the water property data at a particular depth, 
the values obtained depending on the state of the tide at the time of the 
measurements. These fluctuations can make data interpretation very difficult, 
or even lead to mis-interpretation. This could be of considerable significance 
in examining the effects of ocean dumping at an existing dumpsite, or in 
selecting a location for a new dumpsite on the basis of historical data, in 
any area of non-homogeneous water structure. 


A mathematical model of internal tides in fjords would be of great 
utility in evaluating their effect on water property data, especially an 
analytical model of simple application. Dr. Buckley's paper, summarized 
below, set out to determine the feasibility of developing such a model. He 
found several unavoidable complications, including difficulties with suitable 
geometrically-simple representations of the basin shape and with establishing 
suitable boundary conditions. While not completely successful in developing 
an applicable model, his efforts did point out some problem areas, indicated 
the likely region of generation of the internal tides, and substantiated 
the large amplitudes of motion to be expected with this phenomenon in waters 
structured like those of Alberni Inlet. 


Introduction 


The passage of the tide over the sill in a fjord may generate internal 
waves of significant amplitude. Such waves may move water tens of metres 
vertically within the water column and so could lead to erroneous 
interpretation of oceanographic data. Previous studies in Alberni Inlet have 
shown semi-diurnal isotherm excursions of up to 15 m in the inner basin 
(Buckingham, 1978). The present study was designed to provide an analytical 
model of this internal tide and to compare it to existing observations. 


The model, based on an oceanic tidal model of Prinsenberg and Rattray 
(1975), simplified the basin geometry but uses the observed density 
distribution. Because of the geometrical simplifications, the model does 
not exactly reproduce the observed internal tide at all locations in the Inlet, 
but, as will be shown, it does reproduce much of the basic nature of the tide. 


The Model 


The basis of this internal tide model is the standard set of equations 
describing fluid flow. The equations are linearized, inviscid,and two 
dimensional with the effects of the earth's rotation neglected. These 
equations are solved in each of three straight-sided, flat-bottomed regions: 
the outer basin, the inner basin and the harbour. The deep inner and outer 
basins are joined at a shallow sill. There is no obstruction between the 
inner basin and the shallow harbour region. The solutions to the equations 
in the three regions are forced to be continuous at the junctions between 
the regions. 


The fact that the regions are flat-bottomed allows the equations to be 
separated into two parts; one part containing all the horizontal variation 
and the other ‘all the vertical variation. Each of these parts is solved 
independently. The vertical equation is solved by an infinite number of 
functions whose sum represent the total effect of the measured density 
profile. Each of the vertical functions (modes) has a horizontal equation 
associated with it. These horizontal equations are satisfied by simple wave 
solutions whose amplitudes are determined by the basin geometry, boundary 
conditions and region matching conditions. All of the equations for the 
wave amplitudes are solved simultaneously in one large matrix system. The 
matrices are, in principle, infinitely large, but in practice the solutions 
are approximated with a small number of modes. 


Solutions in Alberni Inlet 


Alberni Inlet was parameterized in three regions: an outer basin of 
100 m depth and infinite length extending seaward, an inner basin of 100 m 
depth and 14 km length, and, a harbour region 15 m deep and arbitrarily long. 
A sill extending to 40 m from the surface is located at the junction between 
the inner and outer basins. In both the inner and outer basins, internal 
waves were assumed to travel both towards the head and the mouth of the Inlet; 
in the harbour region, no seaward progressing wave was assumed to exist. 


The density profile used for the test of the model was the average of 
25 profiles taken by Buckingham (1978) in a 25-hour period November 21-22, 
1977. It is shown, along with the first three modes calculated from it, in 
Figure 8. Also shown are the vertical gradients of the first three modes. 
Horizontal currents are proportional to the modes while vertical displacements 
are proportional to the gradients of the modes. The proportionality constants 
are determined from the heights of the observed surface tidal constituents in 
the Inlet and from the tidal prism transports calculated from them. The 
following table shows the amplitudes of internal waves of the first three 
modes at 60 m depth for the two most important tidal contituents. The sum of 
amplitudes for each mode represents the amplitude of internal wave that would 
be observed if only the first three modes were important and is indicative of 
the magnitude of internal wave response. These results show that the observed 
isotherm displacements in Alberni Inlet can certainly be generated by the 
interaction of the surface tide with the bottom topopraphy. 
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Tidal 

Constituents Mode l(m) Mode 2(m) Mode 3(m) Sum (m) 
M» 9.34 -7.44 -12.73 -10.83 
Ky oes 8 0.47 ceo wae, Beret 


The solution to the horizontal equation indicates the source of the 
internal tides. Amplitudes of the horizontal waves for the Mo tide are as 
follows: 


Mode Outer Basin Inner Basin Harbour 
up-inlet down-inlet up-inlet down-inlet up-inlet 

0 1.00 -0.9) 135 +1. 26 29,36 

l 0.00 -1.24 -0.62 -0.63 = One 

3 0.00 O29 0.14 Ole - 0.24 

3 0.00 0.40 0.20 0.20 0.00 


The fact that the up-inlet and down-inlet amplitudes in the inner basin 
are almost equal indicates that almost all the energy is reflected at 
the head of the inlet and therefore the flow of energy from the incoming 
surface (zeroth) mode to the higher modes is caused at the sill. This table 
also shows that the inlet responds more strongly to mode 3 than to mode 2. 
The reason for this result is that the gradient of mode 3 has a zero very 
near sill depth and that of mode 2 does not. 


The horizontal and vertical results may be put together to generate 
a map of internal wave amplitude anywhere inside the model basin, but this 
was not done since there are no data to compare the map with. 


The data that do exist show a tendency for the internal wave amplitude 
to increase with depth, while this model predicts exactly the opposite. A 
region with a sloping bottom has the property of increasing internal wave 
amplitude with depth, and may be a more suitable parameterization of the 
shape of the inner basin. Such a region, however, cannot easily be put into 
the sort of analytic model described here since the vertical and horizontal 
variations do not separate. Explicit inclusion of this shape of basin was 
therefore beyond the scope of the present study. 


Conclusion 


The model has shown that interaction of the surface tide with the sill 
can create internal waves of large amplitude, similar to the observed isotherm 
displacements. It did not predict, however, the correct amplitude variation 
in deeper water. 


The model could be improved by the inclusion of a wedge-shaped region 
as the inner basin and by the calculation of a larger number of modes to 
better approximate the infinite sum solution. A more comprehensive set of 
measurements needs to be made to test the revised model effectively. 


While the model developed here will not predict the characteristics of 
the internal tide at any specific locations in the inner basin of Alberni 
Inlet, it still will be useful in predicting how much energy is taken from 
the surface tide and distributed into the various internal modes under different 
stratification conditions. 


Group Discussion 


In considering inlets with a sloping bathymetry Macdonald asked if 
there was some area where the internal wave would break. Buckley responded 
that such an area likely exists but that the present linear model cannot 
be used to make any prediction about the location of the region. A non-linear 
model would be required for such a prediction. Levings questioned whether 
an internal wave would go through a narrow notch in a manner similar to the 
way it would go through a wide notch. By analogy to light travelling through 
a notch Buckley though that a three-dimensional model would be necessary as 
opposed to the two dimensional used in this study. 
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IX. GROUP DISCUSSION ON FUTURE 
OCEAN DUMPING RESEARCH 


A group discussion of future ocean dumping research trends followed the 
reporting of current projects. 


Studies in Alberni Inlet were summarized and future research projects 
that developed from the Alberni work were noted. Biological data from the 
past few years will be collated and statistically analyzed as part of the 
1979-1980 contracts. This will complete the studies on marine benthos 
communities in Alberni Inlet. 


The chairman noted that there would be some support from the ocean 
dumping budget for the unsolicited proposal submitted to continue a more 
detailed study of the oxygen budget and internal tides of the inner basin of 
Alberni Inlet. This work could complete the need for any further physical 
oceanographic research in the Inlet. 


Chemical research will continue to concentrate on the development of 
meaningful criteria and the establishment of realistic allowable concentrations 
or levels for the various chemical parameters as described in the Ocean 
Dumping Control Act. 


The current status of dumping in Porpoise Harbour and the urgency of 
establishing a more environmentally suitable dumpsite near Prince Rupert were 
discussed. A site-specific examination of the Prince Rupert problem could 
provide information of general utility in the selection of dumpsites. The 
consensus of the group was that research on fish populations and/or 
communities should be supported by Ocean Dumping funds only when such studtles 
were within the framework of studies considering all oceanographic and 
biological features of a study area, and related particularly to dumpsite 
selection or relocation. 


The workshop was adjourned by Mr. Smith. 
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FIGURE 7, COLLECTION STATIONS FOR Molpadia intermedia AT THE 
POINT GREY DUMPSITE. 
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FIGURE 8. PROFILES USED FOR THE MODEL TEST: 
The sigma-t (0+) profile (ain 
Internal mode ( ) and its corresponding 
vertical gradient (----- ). Lon, the firse (bh), 
second (c), and third (d) internal modes. 
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TABLE 3 
SELECTIVE CHEMICAL EXTRACTIONS 


Metal concentrations (ug g7! dry weight) in the different compartments before 
and after the release experiment. Compartments were defined according to the 
selective extraction scheme used by Brannon et al.,19/6. 


COPPER CADMIUM 
VICTORIA HBR. FALSE CREEK VICTORIA HBR. FALSE CREEK 
Before After Before After |__Before __ After _ Before _ After 
IW <0.001* 0.002* IW <0.001* <0.001* 
IE 0.019 0.47 0.015 Lai HES G-02007 0.14 <0.005 0.14 
ER <0.07 O13 <0.05 0.09 |ER 0.068 0.48 <0.05 0.14 
0 133 144 66.7 loz (0. SROs 0.84 1.47 0.50 
MR 2.9 ia | 0.71 <Q.17 |MR <0.068 <0.11 <0%5 <0.08 
RE 29.4 ahs 19.7 40.7 |RE <0.74 <0.95 — <0.58 0.31 
: 165 193 87 B92 Ln eae <1 382 <2.61 Ss doe 
LEAD ARSENIC 
VICTORIA HBR. FALSE CREEK VICTORIA HBR. FALSE CREEK 
Before After Before After |___ Before After _Before _ After 
IW <0.001* <0.001* IW 0.01* <Q.005* 
E 1.03 <253 <0.5 <1) OSes ousore <OV1) <05025 <0.08 | 
ER 4.0.09 0.22 <0.05 <0, 022 1ER <0.34 0.16 <0025 <0.08 | 
0 167 170 38.2 Geen | moto LksG2 La8S 1.00 
MR <6.7 oe <5 1,7 4iMR <0.34 <0.11 <025 <0.08 
RE [232 44.3 29.4 fico Re. ON O2i Awe 0 sSS 0.80 1.88; 
yy 198 222 67.6 62 i Seo <2 0) “S18 <3. 12 
MERCURY 
VICTORIA HBR. FALSE CREEK 
wai BETO FO ATL E a Before. After: KEY 
IW <0.0005* <0.0005* IW Interstitial Water 
E  <0.004 S05 11 <0.003 <0.08 E Exchangeable 
ER <0.034 <Op-id <0.03 <0.08 ER Easily Reduced 
O° «0:034 ee Tt <0.03 <0.08 QO Oxidizable 
MR - AW ealg | <0.03 <0.08 MR Moderately Reducible 
RE {2.14 2.98 0. Sa 0.56 RE Residual 


a on = in 
mie 2221 <3.42 <0.56 <0. 88 Interstitial water units mg L 
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TABLE 4 


SUMMARY OF ANALYTICAL RESULTS FOR AROCHLORS 
OF VICTORIA HARBOUR DREDGE SPOILS 


ANALYSIS MEAN ng g7! RANGE ng g™! 


DRY WEIGHT DRY WEIGHT 
Arochlor 1242 85.8 Ore 210 
Arochlor 1254 205.4 0 - 1950 
Arochlor 1260 S754 0 - 458 
Decachlorobipheny] 11.4 Ora) 337 
Total 360 0 - 3640 
Total (E.P.S. permit 
application 4443-0375) 390 ,84 
TABLE 5 


RESULTS OF STATISTICAL ANALYSIS OF THE PERCENTAGE 
OF VARIATION OF AROCHLORS DUE TO DIFFERENT CAUSES 


Measurement Variation Variation accounted Variation accounted 


| 

| accounted for for between samples for between samples 
| within sub- taken from the taken from 

samples from same barge different barges 

| the same grab 


Sample on the 


barge . 
ne 2 hee ig 


(In %) (in- %) P< (Sin) 
rochlor 1242 44.8 26.8 0225 2625 0.001 
rochlor 1254 40.6 2240 F205 BY PRS 0.001 
rochlor 1260 44.2 26.4 0325 29.4 0.001 
Decachloro- 41.5 45.1 0.001 La4 0.001 
ipheny] 


otal 40.0 a) 0.10 3050 0.001 
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APPENDIX I 


OCEAN DUMPING WORKSHOP ATTENDANCE LIST 


G.R. Smith, Ocean Engineering 1.0.5. (Workshop Chairman) 


C.D. Levings, W.V.L., West Vancouver (Workshop Secretary) 


E.P. Anderson, Dobrocky SEATECH Limited, Victoria 
W.H. Bell, Coastal Zone Oceanography, I.0.S 

Tt. Birtwell’ Habitat Protection, D.<F.0., Vancouver 
W.R. Buckingham, Dobrocky SEATECH Limited, Victoria 
J.R. Buckley, Seakem Oceanography Limited, Sidney 
S.C. Byers, Dobrocky SEATECH Limited, Victoria 

J.A. Crerar, Beak Consultants Limited 

L. Dawson-Burns, MacMillan Bloedel Limited, Vancouver 
R. Deverall, Can Test Limited, Vancouver 

R. Harbo, Habitat Protection, D.F.O., Vancouver 

J. Helbig, Beak Consultants Limited, Richmond 

R.A.W. Hoos, E.P.S., Vancouver 

RAW. Macdonald, Ocean Chemistry, 1.0.5. 

E. McGreer, E.V.S. Consultants Limited, North Vancouver 
D. Morse, Chemex Labs Limited, North Vancouver 

J. Park, Can Test Limited, Vancouver 

S. Steele, D.F.0O., Vancouver 

P. Thomas, Econotech Services Limited, New Westminster 
J.A.J. Thompson, Ocean Chemistry, I.0.S. 

M. Waldichuk, W.V.L. West Vancouver 

C.S. Wong, Ocean Chemistry, I.0.S. 

R. Young, MacMillan Bloedel Limited 


See Appendix III for details of the research institution abbrevlattons. 
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APPENDIX II 


1978-79 CONTRACTS 


Blolopy 
| Ref. DSS File No. O08SS.  KF833-7-1210 $27,000.00 
- Ref. DSS File No. O8SB. KF833-8-0248 $29,010.00 


Marine biological sampling and observations of temporal changes 
in benthic communities and benthic respiration at a dumpsite in 
Port Alberni. 


Scientific Authority: C. D. Levings, W.V.L. 


Contractors: lReak Consultants Ltd. 


2Dobrocky SEATECH Ltd. 


Ref. DSS File No. O8SB. KF833-8-0248 
Amendment a oles U0 


(a) Review Iiterature on the reproductive ecology of selected 
invertebrates from a dumpsitte in Alberni Inlet. 


(b) Verify identification of invertebrates performed by others 
working at Alberni Inlet and compare identification with current 


studies. 


Scientific: Authority: C. D., bevings, W.V.L. 


Contractor: Dobrocky SEATECH Ltd. 
Chemistry 
Ref. DSS File No. 07SB. KF833-8-0590 $ 6,300.00 


Analysis of sediment samples for polychlorinated biphenyls. 
Scientific Authority: R. W. Macdonald, I.0.S. 
Contractor: Seakem Oceanography Ltd. 

Ret, 205s PLIezNO. O7SB., (KF835-38-0591 > 7,390.00 


Determination of the lead methylating capacity of Pacific coastal 
marine sediments. 


Scientific Authority: J. A. J. Thompson, I.0.S. 


Contractor: Beak Consultants Ltd. 


C. 


ip ES 


Ref, DSS File No. O7SB. KFS33-d=0592—1 SLO,421. 05 


Development of a meaningful criteria for ocean disposal of dredged 
or sedimentary material - Part l. 


Scientific Authority: R. W. Macdonald, [.0.S. 
Contractor: Can Test Ltd. 
Ref. DSS File No. O7SB. KF833-8-0592~2 S10e241G25 


Development of a meaningful criteria for ocean disposal of dredged 
or sedimentary! material = Part II. 


screntifie*Authorrry sk. iWy Macdonald; t.O°S. 
Contractor: Seakem Oceanography Ltd. 
Ref. DSS File No. O8SB. KF833-8-1577 > pie U 


Analyses of samples of marine holothurian, Molpadta intermedia 
for the determination of copper, zinc and cadmium. 


ScheneLerc. ANENOLICyY: od. As J, LNOMpsONsa | .O.0. 


Contractor: Can Test Ltd. 
Physics 
Ref. DSS File No. O7SB. KF833-8-0594 S$ 7,980.60 


Examination of the variability of upwelling on the West Coast 
of Vancouver Island and its relationship to the flushing of 
Alberni Inlet. 

sereltiiteAurtiority: W. Ns Belli T.0%S. 

Contractor: Beak Consultants Ltd. 


Ref. DSS File No. O7SB. KF833-8-0593 > 4,402.00 


Examination of evidence for an internal tide in Alberni 
Inlet. 


Scientitic Authority: “Wa N. Bello Eo OeS. 


Contractor: Seakem Oceanography Ltd. 


Bie ae 


Administrat fon 

Ret. DSS File No. O7SB.  KF833-8-0459 S 8,435.00 
Support for the West Coast dumping program. 
Scientific Authority: R. O. Brinkhurst, I.0.S. 
Contractor: Dobrocky SEATECH Ltd. 

For H.Q. Office 52,900.00 


Annotated bibliography of Canadian papers and reports bearing on 
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ABSTRACT 


Infrared satellite imagery of the waters off the Pacific coast of Canada 
for the five-year period, July 1974 through July 1979, has been examined in 
order to delineate sea surface temperature patterns in the area. A series of 
charts depicting oceanic thermal fronts is compiled. The source data consist 
of enhanced VHRR, NOAA satellite infrared imagery in photographic-print format, 
obtained mostly from the National Satellite Data Services Branch of the 
National Climatic Center (NOAA) and some from its field station (San Francisco, 
California) and also from the Arctic Weather Centre of the Atmospheric Environ- 
ment Service of Canada (Edmonton, Alberta). Some are based on original data 
acquired from the San Francisco field station and the Arctic Weather Centre. 
With regularly-obtained imagery in photographic-print format only two shades 
of the grey-tone variations of the sea surface can generally be delineated 
successfully by an optical device. However, with suitably enhanced imagery 
made from original data, at least six shades of these variations can be dif- 
ferentiated. Unless enhancement is made consistent from image to image it is 
possible to misinterpret locations of weak fronts. In order to track fronts 
successfully from day to day it is desirable to use original satellite data. 
Comments are made to explain the occurrence of certain water types at various 
locations. 


Introduction 


Sea surface temperatures obtained by the infrared radiation sensors aboard 
satellites have been applied, since 1968, to depict a number of oceanographic 
features and processes. The scientific applications range from studies of up- 
welling (e.g. La Violette and Chabot, 1968; La Violette, 1974; Stumpf, 1975) to 
those of ocean circulation (Vukovich, 1974); eddies (e.g. Gotthardt and Potocsky, 
1974) and oceanic fronts (Wyrtki, 1977). Techniques for processing satellite 
data have advanced sufficiently to a stage that now daily or weekly charts of 
distribution of sea surface temperatures over wide oceanic areas are compiled 
routinely and used operationally by countries such as Japan (Hanzawa, 1979) and 
the U.S.A. (Pichel et al, 1979). Moreover, detailed charts of fishing areas are 
produced regularly by France (Noel, 1979), Japan (Hanzawa, 1979) and the U.S.A. 
(Jurick, 1977) and made available to commercial fishermen as aids to their fish- 
catching effectiveness. 


Over the past five years (1974-1979) well over 100 infrared satellite 
imagés of the waters off the Pacific coast of Canada and northern U.S.A., for 
relatively cloud-free conditions have been acquired. Inspection of these re- 
vealed that less than one half of the total number possessed sufficient defi- 
nition of sea surface thermal features necessary to map oceanic thermal fronts 
in the area. The present report is primarily a compilation of charts depicting 
major oceanic fronts off the Pacific coast of Canada (Figure 1) based on these 
images. It also contains some discussion on data quality and on oceanic events 
deduced from the examination of the imagery. 


Data 


Satellite data obtained by the very high resolution radiometers (VHRR) of 
the spacecrafts, NOAA-3, NOAA-4, NOAA-5 and TIROS-N, were utilized. The visual 
and infrared imagery based upon these data are reported to have a spatial re- 
solution of one kilometre at subpoint. The full sensitivity of the infrared 
sensors aboard these satellites is + 0.5°C (Koffler, 1976). The majority of 
the data were procured from the National Satellite Data Services Branch of the 
National Climatic Center, NOAA (Washington, D.C.) while the rest were obtained 
from its San Francisco Field Service Station (Redwood, California). A limited 
number was obtained from the Arctic Weather Centre of the Atmospheric Environ- 
ment Service (Edmonton, Alberta). The imagery are available in 25 cm x 25 cm 
size photographic-print format. Examples of these, in reduced size, are shown 
in Figures 2 and 3. The imagery is essentially photographic prints produced 
by a computer-interfaced Muirhead-type facsimile recorder. The visible data 
(Figure 3a) are displayed over the full range of the reflective surface bright- 
ness, the white and black representing respectively the light and dark colour 
of the surface. The infrared data (Figures 2 and 3b) cover a large range of 
infrared radiation temperatures. Because of this, neither the sensitivity of 
the photographic paper nor the eyes will be able to perceive small changes of 
temperature, say, a degree. In order to overcome this deficiency usually a 
narrower range is selected to give an "enhanced" image covering, for example, 
the range 0 to 20 C in grey-tone levels. Light tones represent lower temp- 
eratures; darker tones represent higher temperatures. The white and black on 
these prints now represents, respectively, the temperatures less than and 
greater than those assigned as limits. In the imagery obtained routinely 
from the Climatic Data Center, even the enhanced data are not all suitable for 
a local area on account of the choice of larger range to cover greater oceanic 


areas. Generally, imagery obtained from the field station contained better 
definition of the sea surface features than did that from the Climatic Data 
Center - a reflection of the difference in enhancement made at the two data 
sources. 


A limited amount of original satellite information which had been digit- 
ized and recorded on magnetic tapes at the San Francisco Field Service Station 
and the Arctic Weather Centre was obtained and the data processed locally. A 
few examples of these are included in this report to show the superiority of 
these products over regularly-obtained photographic prints. 


Data Analysis 


The visual and infrared imagery were viewed in pairs (in the case of day- 
time observation) to determine the extent of cloud-free oceanic areas. Only 
such areas were considered for further examination. Cloud-free areas for 
nighttime imagery were determined from interpolation of two successive daytime 
observations, one before and the other after the observation in question. Data 
that presented difficulty with this interpolation on account of possible fog or 
low cloud formation during the night were excluded from the analysis. Areas 
that did not possess distinguishable surface features as a result of poor data 
quality, such as inadequate imagery enhancement and/or "noisy'' scan-lines, etc., 
also were excluded from analysis. 


An optical device, the Zoom Transfer Scope (Bausch and Lomb), was utilized 
to transfer the thermal surface features, such as an oceanic front, in the 
imagery, to a base chart. This device can readily reduce or enlarge imagery 
to an appropriate scale and its anamorphic system can compensate for geometric 
distortion in the imagery due to relief, earth curvature, print shrinkage, etc., 
thereby allowing a feature on a print to be superimposed onto a chart or map. 
When control points on the imagery are clearly distinguishable, such as small 
islands, sharp points or capes of coastline, positions of sea surface features 
in their general vicinity can be transferred to a base chart with an estimated 
accuracy of a few kilometres. In areas far removed from such control points 
the accuracy is expected to drop considerably, by up to an order of magnitude. 


In some prints the area under study was located at the edge or corner of 
the prints where large distortion of the earth's surface was present. In such 
a case the Zoom Transfer Scope is incapable of correcting the spatial distor- 
tion. For these a transparent grid overlay supplied by NOAA and containing 
latitude and longitude lines was used to obtain positions of sea surface 
features. The error in the determination of the positions using this method 
may be as large as 30 km (Stevenson and Kirkham, 1977). 


The original data sets recorded on magnetic tape were processed with the 
aid of a Hewlett-Packard 2100S minicomputer and Versatec 1100A printer-plotter. 
The selected parts of the digitized radiance information were printed as grey- 
tone variations on a printer-plotter. The image was geometrically corrected 
for the strong side-to-side distortion and several adjacent pixels (picture 
elements) were averaged to reduce effects of noise and banding. Averages 
of 2 x 4 and 7 x 8 pixels were made, the latter providing information averaged 
over an area of 6 km x 7 km. Oceanographic data collected approximately at the 
time of satellite observations were used to "field-calibrate" the latter. An 
example of the data used for calibration is shown in Figure 4 (Tabata and Gower, 


1979). It is evident from these data that the satellite data are capable of 
yielding sea surface temperatures to within 0.5°C, a figure which is better 

than that of the routinely-observed satellite temperatures of + 2°C (Brown, 

1975). 


Results 


In Figure 5-1 through 5-52 are shown illustrations depicting areas of 
warm and cool waters and their associated fronts off the Pacific coast of 
Canada and northern U.S.A. They were interpreted from infrared satellite 
imagery available in photographic-print format. 


As is evident from these illustrations, during spring through autumn, 
cooler waters were present adjacent to Vancouver Island at three locations: at 
the approaches to and in Juan de Fuca Strait, at the approaches to and in 
Queen Charlotte Strait and in the vicinity of Scott Islands (e.g. 28 April 
1976 = Fig. 5-293) 31' May 1978.= Fig..5=-44; 30 July 1974--— Fig. 5-13 and#20 
October 1977 - Fig. 5-42). The few sets of data for winter show that during 
this season the waters at these locations were warmer than in their immediate 
vicinity (e.g. 5 January 1979 - Fig. 5-47) except on 31 December 1977 (Fig. 
5-43) when it was cooler. The contrast between the cool and warm water is 
most striking during summer and autumn and least in winter and spring. 


The waters present along the west coast of Vancouver Island were usually 
cooler than those lying immediately offshore (e.g. 5 January 1979 - Fig. 5-47; 
4 August 1974 - Fig. 5-3 and 20 October 1977 - Fig. 5-42). However, from time 
to time isolated areas of warm water occurred adjacent to the coast (e.g. 28 
April 1976 - Fig. 5-29; 31 May 1978 - Fig. 5-44; and 9 September 1975 - Fig. 
5-20). 


Generally, the warm water that lay offshore occurred in a band of few 
hundred kilometres breadth at latitude 49°N (e.g. 10 September 1976 - Fig. 
5-20). In one example, however, the band was as narrow as 50 km (31 December 
1977 - Fig. 5-43). Landward of this band of warm water was a well-defined 
tongue of cool water with its source along the west coast of Vancouver Island 
and its leading edge directed southward. This tongue was most easily recog- 
nized during summer and autumn (e.g. 14 July 1979 - Fig. 5-52 and 9 September 
1975 - Fig. 5-18). 


One of the more interesting features of the sea surface temperature 
pattern off the coast, especially off Vancouver Island, is the presence of a 
number of well-defined plumes or tongues of cool water (in most cases) and of 
warm water with their axes lying perpendicular to the coast. They occurred 
during all seasons (e.g. 5 January 1979 - Fig. 5-48; 1 June 1978 - Fig. 5-45; 
4 August 1974 - Fig. 5-3; and 9 September 1975 - Fig. 5-18). The separation 
between the axes of the tongues was typically 100 to 150 km. Some of the 
tongues of cool water extended seaward, intercepting Line P at Stations 3 and 
4. 


The cool water of Juan de Fuca Strait often spread northwestward until it 
finally lost its identity as it merged with the other cool water lying off 
Amphitrite Point (e.g. 26 August 1974 - Fig. 5-5 and 18 October 1976 - Fig. 
5-32). Occasionally one branch of this water spread to northwestward while 
another branch backed and spread southward upon leaving the Strait (e.g. 5 


January 1979 - Fig. 5-47 and 1 June 1978 - Fig. 5-45). At other times this 
northwestward "flow" continued to the vicinity off Amphitrite Point, then 

backed and "flowed'"' southward, intercepting Line P in the vicinity of Station 

2 and occasionally Station 1 (e.g. 23 August 1975 - Fig. 5-16; 21 September 1975 
- Fig. 5-26: and 20 October 1977 - Fig. 5-42). 


The cool water lying off the mouth of Columbia River in winter (6 
February 1976 - Fig. 5-27) and warm water there in summer (31 May and 1 June 
1978 - Figs. 5-44 and 5-45; and 14 July 1979 - Fig. 5-52) is likely to be 
associated with the Columbia River plume. 


The Queen Charlotte Sound - Hecate Strait area is usuallv cloudier than the 
southern areas and consequently there are less data available there. From the 
few sets of data available it appears that a large tongue of warm water with 
its source in the offshore region further south, intruded the area (e.g. 9 
September 1975 — Fig. 5-18; 20 October 1977 - Fig. 5-42 and 14 July 1979 - 

Fig. 5-52). Offshore of McInnes Island the water was frequently cooler than 
the adjacent waters during all seasons (e.g. 5 January 1979 - Fig. 5-47; 1 
June 1978 - Fig. 5-45; 29 August 1974 - Fig. 5-7 and 20 October 1977. - Fig. 
342 x 


An example of the distribution of sea surface temperatures of the study 
area produced from the original data in digital form is shown in Figure 6. The 
series of numerals in Fig. 6a represent the grey-scale levels of the infra- 
red radiation from the earth's surface. By making use of the calibration 
curve (Fig. 3), these levels in turn can be converted into equivalent sea sur- 
face temperature levels from which appropriate isotherms can be drawn. Fig. 
6b shows the detailed sea surface temperature distribution for 9 September 
1975 as derived from the data indicated in Fig. 6a. Another example from 10 
September 1975 is shown in Fig. 7. Despite the fact that the above pair of 
observations was made from two separate spacecraft, the former by NOAA-3 and 
the latter by NOAA-4, the temperatures for the two consecutive days are re- 
markably similar. 


Discussion 


Only one half of the total of more than one hundred infrared satellite 
imagery in photographic-print format on hand possessed sufficient definition 
of sea surface temperature patterns required for mapping oceanic fronts. Of 
this. some still presented difficulty with imagery interpretation due to fuzzy 
definition and "noisy" scan lines. The fact that some of the imagery con- 
tained featureless sea surface does not necessarily imply that features were 
absent. Rather, it would appear that this could be the result of the degree 
of enhancement made on each data set. Imagery obtained from the Climatic 
Data Center seems generally to be enhanced over a relatively large temperature 
range. This type of enhancement is not desirable for certain local areas 
where more effective enhancement can be made over a narrower temperature range. 
In the results presented in Fig. 5 it is noted that most of the fronts shown 
for successive days are generally consistent from day to day. On the other 
hand, there are some that portrayed appreciable variations such as during 26 - 
28 November 1976 (Fig. 5-34 to 5-39). Some of the variations might be real, 
due to day to day changes in surface conditions, but others are believed to 
be due to the degree of enhancement made on each set which would have the 
effect of shifting the position of the fronts. 


It was initially thought that it should be possible to delineate at least 
four shades of the grey-tone variations in the photographic-print imagery 
employing optical devices. However, this proved to be more difficult than 
was originally assumed. For most, only two shades could be delineated proper- 
ly. In contrast to this, with imagery that had undergone enhancement to suit 
local conditions it was possible to differentiate at least six shades of the 
grey-tone variations, as has been possible for the data of 14 July 1979 (Fig. 
5-52). 


The few examples of detailed sea surface temperatures obtained from the 
analysis of raw data demonstrate the effectiveness of using original data to 
map local sea surface temperature patterns. The remarkable similarity between 
the patterns obtained from NOAA-3 during its nighttime pass of 9 September 
1975 (Fig. 6b) and that derived from NOAA-4 during its daytime pass of 10 
September 1975 (Fig. 7) is an indication of its effectiveness over patterns 
derived from routinely-obtained imagery from the data centre. 


The imagery revealed a variety of features in the distribution of sea 
surface temperatures off the coast upon which some comments can be made. 


The occurrence of cool water along the Pacific coast of North America 
during the summer months is usually attributed to the effect of upwelling 
(Sverdrup et al, 1942). The northerly winds that prevail along the Pacific 
coast during spring through autumn result in the divergence of surface water 
adjacent to the coast as a result of the offshore Ekman transport. In order 
to compensate for the surface offshore drift, cold and saline water is up- 
welled from greater depth. Off Oregon upwelling can occur as soon as early 
spring but off the Pacific coast of Canada it does not appear until late 
spring (Bakun, 1973). Studies conducted over the continental shelf of Oregon 
indicate that the region of upwelling is restricted to a very narrow band of 
approximately 10 km from the edge of the coast (Halpern, 1974; Moores” etal ; 
1976). For the mean oceanographic conditions over the continental shelf off 
Vancouver Island the Rossby baroclinic radius of deformation is approximately 
40 km; consequently, the region of upwelling here is likely to extend further 
offshore (Yoshida, 1955). 


The frequent occurrence of a band of cool water adjacent to Vancouver 
Island and the Washington-Oregon coast, as determined from satellite imagery, 
can in many cases, be explained in terms of upwelling. For example, the 
persistent occurrence of cool water during 30 July - 5 August 1974 (Figs. 5-1 
to 5-4) is almost certainly due to the effect of upwelling as the upwelling 
indices for the area were favourable during the one month between late July 
and late August (Bakun, in preparation). During other periods upwelling con- 
ditions were not as persistent for extended duration. Such were conditions 
during late August and September 1974. The rather "spotty" occurrence of 
cool water adjacent to the coast during this period (26-29 August 1974 - 

Figs. 5-5 to 5-7; 10-13 September 1974 - Figs. 5-8 to 5-11; and 19-20 September 
1974 - Figs. 5-12 to 5-14) is probably a reflection of irregular upwelling 
along the coast. Some of the imagery showed that, even when upwelling con- 
dition was favourable, as indicated by the indices for latitude 48°N, longitude 
125°W and latitude 51°N, longitude 131°W, the band of cool water did not extend 
continuously from Vancouver Island to the Oregon coast (28 July 1975 - Fig. 
3-15; 23 August 1975 - Fig. 5-16; 8-12 September 1975 - Figs. 5-17 to 5-22; 


and 18-21 September 1975 - Figs. 5-23 to 5-26). The reason for this discrep- 
ancy may be attributed either to the indices not being representative for the 
local area or else to certain areas having experienced local heating. 


The occurrence of cool water along the coast during the winter (31 
December 1977 - Fig. 5-43; 5-7 January 1979 - Figs. 5-47 to 5-49 and 21 
January 1979 - Fig. 5-50) is unlikely to be due to the upwelling effect. For 
the case of 31 December 1977, it may be due to cooling from excessive heat 
loss to the atmosphere that was associated with the local atmospheric circu- 
lation at that time. Starting on 29 December very cold continental air was 
blowing offshore for more than a week (Atmospheric Environment Service, 
private communication), In such a situation it can be assumed that the largest 
cooling over water occurs closest to the land. It is also possible that mist, 
which might give misleading sea surface temperature in the infrared imagery, 
could have formed over the water, but none of the coastal weather stations 
reported it. For the other two periods in January 1979, cooling due to heat 
loss is also suspected as upwelling condition was not favourable during these 
periods (Bakun, in preparation). 


The cool water associated with the tongues has lower salinity than 
does the adjacent water (Tabata and Gower, 1979) and therefore can be considered 
as the water that originated near the coast. This would imply that the cool, 
brackish coastal waters are advected seaward as relatively narrow jets rather 
than in a broader diffused flow. 


The cool waters that are present during spring through autumn and warm 
waters there during winter at the approaches to and in Juan de Fuca Strait 
and Queen Charlotte Strait and in the vicinity of Scott Islands can be attri- 
buted to tidally-induced vertical mixing. These areas are characterized by 
the presence of strong tidal currents with their speed exceeding 1's m/sec at 
certain parts within them (Environment Canada, 1976). Even in areas where 
tidal currents are not as strong, there appears to be recognizable tidal in- 
fluence on sea surface temperatures observed in inland seaways due to tidal 
mixing (Herlinveaux, 1957). During late spring through early autumn the 
coastal waters are generally warmer at the surface than immediately below and 
vice-versa in winter (Herlinveaux and Giovando, 1969). The consequence of 
tidal mixing would then be to reduce the surface temperature during spring 
through autumn and to raise it during winter. It is to be noted that wind 
mixing could achieve the same results also. 


Vertical mixing capable of overcoming stratification of coastal waters 
need not be restricted only to waters adjacent to land. Indeed, such mixing 
can be effective in more open waters of the continental shelf if the current 
speeds there are sufficiently large. In this regard, Simpson and Hunter (1974) 
have proposed that the transition between stratified and mixed regimes can be 
represented by a mixing parameter based on the ratio of potential energy 
required for complete vertical mixing to the rate of dissipation of energy of 
the tidal currents. Where the buoyancy term remains relatively constant the 
parameter simplifies to H/U3, where H is the water depth and U is the ampli- 
tude of the tidal current. This, of course, assumes that other currents are 
of minor importance. For the continental seas off the United Kingdom, com- 
plete vertical begs: was found to occur when the value of this parameter was 
less than 50-100 sec?/m?. 


We can apply this criterion if we assume comparable tidal friction and 
water column stability for the waters off the Canadian coast. Over the conti- 
nental shelf off Vancouver Island, the amplitude of the tidal current is about 
50 cm/sec, based on only one series of observations off Tofino (Huyer, et al, 
1976). From a very limited number of observations in Queen Charlotte Sound- 
Hecate Strait (Canadian Joint Committee on Oceanography, 1955) it would appear 
that the amplitude of the tidal currents herein is of the same order of magni- 
tude as off Vancouver Island. Given an amplitude of 50 cm/sec, only areas with 
depth less than approximately 10 m would experience complete mixing upon apply- 
ing the criterion. Since very little area of the shelf off Vancouver Island 
and Queen Charlotte Sound-Hecate Strait is as shallow as 10 m, tidal currents 
with the above amplitude would be ineffective in destroying the stratification 
necessary to lower the surface temperature in summer and raise it in winter. 
However, if the tidal amplitude had been underestimated and if we can assume 
that the amplitude can be 50% greater than quoted above, areas with depth less 
than 30-40 m should experience complete mixing, even in the absence of other 
currents, or of wind or convective mixing. Such an area is present along the 
western half of northern Hecate Strait and a small area (15 km x 20 km) loca- 
ted about 80 km south of McInnes Island. The occasional occurrence of cool 
water along the western side of Hecate Strait (e.g. 19-20 September 1974 - 
Figs. 5-12 to 5-14 and 9 September 1975 - Fig. 5-18) could have been due to 
tidally-induced vertical mixing. No evidence of cool water in summer has been 
found for the area south of McInnes Island. Probably the tidal current there 
is not sufficiently strong to overcome stratification. 


The presence of cool water at the approaches to Juan de Fuca Strait can be 
explained on the basis of other than wind-induced upwelling and tidal mixing men- 
tioned earlier. It can also be a consequence of local upwelling associated with 
the behaviour of the seaward flow of Juan de Fuca water (Tully, 1942). This 
follows from the assumption that the seaward flow is in the form of a "wake 
stream" or a jet. This stream, in northern hemisphere, should veer to the 
right and a counter current should occur along the left side of the wake 
stream (Rossby, 1937-38). Such a current system could result in a formation 
of a zone of upwelling between these two currents. 


In order to ascertain which processes are responsible for the occurrence 
of cool water at the approaches, other supporting information such as salinity 
would be required. As has been shown recently (Tabata and Gower, 1979), only 
with salinity data was it possible to interpret the cool water there on 9-10 
September 1975 (Figs. 5-19 and 5-20) as being due to local upwelling. 


The Columbia River water, before entering the ocean, is warmer than the 
ocean water at the river mouth in the summer and cooler during the winter, by 
a few Celsius degrees in both cases (National Ocean Survey, 1970). According- 
ly, during these seasons at least, the Columbia River plume should be clearly 
distinguishable from the oceanic waters. Yet, of the 27 images that included 
the area off the mouth of Columbia River only a few indicated its presence 
with any certainty. The reason for its absence in the imagery for the summer 
months is not clear. One possible explanation is that the warm river water 
upon entering the ocean is rapidly cooled by mixing with the cool, upwelled 
waters along the coast. It is also possible that the degree of enhancement 
was insufficient to delineate the plume from the surrounding waters. 


During the winter the Columbia River plume appears to flow northward in a 


narrow band of few tens of kilometres off the coast (Budinger et al, 1964). 
Accordingly, even if the plume water were cooler than the surrounding waters 
it will be difficult to delineate the plume from the other coastal waters that 
might have been cooled by other means, such as heat loss to the atmosphere. 
The band of cool water present along the Washington coast during 31 December 
1977 (Fig. 5-43) and 5 January 1979 (Fig. 5-48), at least for the portion 
close to the river mouth, might have been the plume but from temperature data 
alone it will be difficult to ascertain this. Salinity data should be able 

to resolve this question. 


The occurrence of a number of tongues of cool or warm water with their 
axes frequently oriented perpendicular to the coast and spaced regularly (e.g. 
Fig. 5-18) suggests the presence of wave-like phenomena. Mysak (1977) has 
examined the above feature and proposed that the observed pattern may be due 
to the baroclinic instability of the California Undercurrent. Detailed analy- 
sis of satellite data obtained during the two days in early January 1979 has 
provided further evidence that the cellular pattern evident in the satellite 
imagery is consistent with the model of baroclinically unstable waves (Emery 
and Mysak, 1979). Whether or not such pattern can definitely be attributed to 
this effect would require additional supporting data. 


Conclusions 


Only one half of a set of satellite imagery obtained from the Climatic 
Data Center could be interpreted adequately. The lack of features of the sea 
surface does not necessarily imply that the sea surface was featureless. In- 
stead, the absence of features is almost certainly to be due to the degree of 
enhancement applied to each set of data. Unless the degree of enhancement is 
consistent from image to image it is even possible to misinterpret the sequence 
of features; for example, difference in enhancement can shift positions of 
some of the less well-defined fronts. With the regularly-obtained imagery 
from the Climatic Data Center generally only two shades of the grey-scale 
level could be delineated with a conventional optical device. On the other 
hand, it is possible to differentiate at least six shades of the grey-scale 
level in a suitably-enhanced image. Due to difficulty associated with en- 
hancement, daily monitoring of movement of fronts from successive imagery 
could be misleading. For any detailed analysis of sea surface temperature 
pattern, especially for the tracking of the movement of fronts and eddies in 
local areas,'such as the continental shelf, analysis of original data is 
preferable as the enhancement can be made consistent for any sequence of 
imagery or it can be varied to suit the task on hand. For coastal waters the 
use of satellite imagery alone to interpret oceanographic processes can be 
misleading as there may be a number of processes that may change water tem- 
perature. Other properties of water, such as salinity would be required to 
validate the interpretation. 


The oceanic region off the Pacific coast of Canada is frequently cloud 
covered. For about 90% of the time either clouds or fog are present. As a 
result useful infrared satellite coverage over the area is limited. However, 
despite this, from the examination of satellite imagery obtained over the 
past five years it has been possible to get some better insight into the 
nature of circulation off the coast. With improved access to the original 
satellite data it should be possible to use the data effectively, in "real 
time", not only as aids to the study of oceanographic processes and events 


and even in the planning of oceanographic surveys, but also as a guide to 
commercial fishermen searching fishing grounds and to those involved in 
coastal pollution study and pollutant cleanup. 
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(b) Isotherms (°C) constructed from data shown in Fig. 6a. 
The dashed lines over the ocean represent the cruise 
track of CFAV Endeavour (19 August - 10 September 1975). 


Sea surface temperatures (°C) derived from original infrared 
satellite data of 1709 GMT, 10 September 1975 (these data 
were used by NOAA to produce the imagery shown in Fig. 3). 
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Figure l. Chart showing general area of study. The weathership obser- 
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indicated. 
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Figure 2. Very high resolution radiometer (VHRR) NOAA-3 infrared 
imagery of western Canada and U.S.A. and of the waters 
off their Pacific coast, based on data taken during the 
nighttime overpass at 0339 GMT, 9 September 1975. The 
imagery has been enhanced to accentuate sea surface thermal 
features. The lighter and darker tones of the ocean are 
associated with cooler and warmer waters, respectively. 
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Very high resolution radiometer (VHRR) NOAA-4 satellite 
imagery of western Canada and U.S.A. and of the waters 
off their Pacific coast, based on data taken during the 
daytime overpass at 1709, GMT, 10 September 1975: 


(a) Visual imagery, showing among other things, the cloudy 
and cloudfree areas. The white and black represents 


respectively the light and dark colour of the earth's 
surface. 


Figure 3. 
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Very high resolution radiometer (VHRR) NOAA-4 satellite 
imagery of western Canada and U.S.A. and of the waters 
off their Pacific coast, based on data taken during the 
daytime overpass at 1709, GMT, 10 September 1975: 


(b) Infrared imagery (enhanced). The lighter and darker 
areas of the ocean are associated with cooler and 
warmer waters, respectively. 
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shipborne and shore-station sea surface temperature. The 
satellite data were taken at 0339 GMT, 9 September 1975; 
shipborne and shore-station data were collected within one 
day of the satellite observation. The location of shore- 
stations (Race Rocks, Amphitrite Point, Kains Island and 
Pine Island) is shown in Fig. 1. Shipborne data were col- 
lected near latitude 48°N and within 500 km of the coast. 


Figure 5. 


ve 


Chart of the ocean off the Pacific coast of Canada and north- 
ern U.S.A. showing areas of relatively warm and cool waters 
and the fronts associated with them, as deduced from satellite 
imagery. The continuous and dashed lines delineating the 
warm and cool water represent strong and weak fronts, respect-— 
ively. The blank space can be considered either as an area 

of no data coverage or where the definition of features was 
poor. Each illustration has annotated data: time (hours, 
minutes, seconds, GMT), date, orbit number and spacecraft 
identification. 
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Figure 6. Sea surface temperatures derived from original infrared 
satellite data of 0330 GMT, 9 September 1975 (these data 
were used by NOAA to produce the imagery shown in Figure 2). 


(a) Grey-scale digitized imagery. The numerals represent 
the averages of 7 x 8 pixels (equivalent to average 
temperatures of an area 6 km x 7 km), the lower and 
higher numerals representing the lower and higher 
temperatures, respectively. The area covered in this 
illustration is shown as the small bordered portion in 
Big. 2; 
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Sea surface temperatures derived from original infrared 


satellite data of 0330 GMT, 9 September 1975 (these data 
were used by NOAA to produce the imagery shown in Figure 2). 


Figure 6, 


(b) Isotherms (°C) constructed from data shown in Fig. 6a. 
The dashed lines over the ocean represent the cruise 
track of CFAV Endeavour (19 August - 10 September 1975). 
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Figure 7. 
satellite data of 1709 GMT, 10 September 1975 (these data 
were used by NOAA to produce the imagery shown in Fig. 3). 
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ABSTRACT 


A static single-point mooring model is developed for use in the 
analysis of moored instrument array geometry and for the interpretation of 
collected data. The model incorporates the steady three-dimensional 
hydrodynamic forces exerted by a space-varying velocity field acting on a 
moored system comprised of a subsurface buoy, a stretchable cable and 
various instrument packages. The system equations are solved by piecewise 
integration down the cable. The computer program provides a comprehensive 
listing of coordinates, angles, forces, etc. 
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Introduction 


The Coastal Zone Oceanography Group at the Institute of Ocean 
Sciences, located at Sidney, B.C., makes frequent use of moored instrument 
arrays for the measurement of ocean currents and seawater properties. The 
moorings are usually single-point taut-line types, often supported by 
subsurface buoys. Some time ago, it was found expedient to develop a static 
two-dimensional model to aid in the design of the moored systems and in the 
analysis of the data obtained therefrom. This development was described in 
Pacific Marine Science Report 77-12 (Bell, 1977), entitled Static Analysts 
of Stngle-Point Moorings. The two-dimensional model has proved to be a 
valuable design tool, but has sometimes been inadequate for problems of the 
“hindcasting" type because the data are, very often, three-dimensional in 
nature. Therefore, it was decided to revise the model to provide a three- 
dimensional capability, as reported herein. 


Certain material included in PMSR 77-12, pertaining to the utility 
of models, mooring cable elasticity and hydrodynamic drag, is not repeated. 
We suggest that a perusal of that report might be a worthwhile preliminary 
to a study of the present one. Although not a prerequisite for proper 
application of the three-dimensional model given here, it could provide some 
additional insight into the problem. The capability for using either a 
surface or subsurface buoy, which was incorporated into the computer program 
for the two-dimensional mooring model, has not been retained in the new 
model. In retrospect, it appeared simpler to write and apply separate 
routines for the two cases. The three-dimensional model does include such 
details as elastic cable, varying cable types along the length of the mooring 
line, instrument or buoyancy packages attachable along the line and a space- 
varying velocity profile. 


Theoretical Analysis Of Ehe Buoy-Cable system 


A. The Buoy 


No particular buoy geometry has been assumed in the development of 
the model. In practice, the shape of the buoy may result in the generation 
of lift forces, or moments of forces, and no account of such complications 
is taken here. The simplest case only is considered, the basic assumptions 
being: 


1. The line of action of all forces is through the centre 
of gravity of the buoy. 


Zz. ‘There are no cross-fhow forces, i.e.fAano lift or kiging 
forces. 


3. There is no change in drag coefficient with a change in 
buoy aspect. 


The justification for retaining as much simplicity as possible is that the 
mooring system configuration is primarily determined by weight and buoyancy 
forces and by cable drag, buoy drag usually being a relatively unimportant 
feature. 


To begin, flow velocity components U, V and W are specified ina 
Cartesian coordinate system (called here the principal reference frame) 
along the directions of the X, Y and Z axes, respectively, where Z is 
vertical and positive upwards. The components (Fig. 1) have a resultant 
velocity vector R whose orientation is defined by the angle $ in the vertical 
plane and y in the horizontal plane, and whose magnitude is given by: 


ne (1) 
Therefore: 
2 a 
R = mh F oe ma rr eee f — ae = i (2) 
sin ¢ sin y sin ¢ cos y cos sin 9 
and 
age ashe AA 
d = cos | 2 | (3) 
y = tan?{ 2 | (4) 


The angle ¢ is here called the current tilt and y is the current direction. 


The buoy is considered as a free body, with the various force 
components acting through its centre of gravity. It is assumed that the 
buoy drag coefficient C is independent of the aspect presented to the 
current, as stated above. The cross-sectional area of the buoy is Ap: The 
drag force is given by: 


Oy aoe . A, R|R| (5) 
with components: 

Day = D, sin > sin y 

Day = Da sin ¢ cos y (6) 

Daz = D, cos 


The weight force W, and buoyancy force Ba have components in the Z-direction 


B 
only, so that the gravitational force components are: 
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(7) 


Tension components, due to the reaction of the mooring cable at the point of 
attachment, must also be taken into account. The orientation of the cable 
(Fig. 1) is defined by the angles 0 and wy, where 9 lies in the horizontal 
plane and is here called the cable azimuth, and w~ is in the vertical plane 
and is called the cable angle. The tension components are: 


Ty = -T sin ) sin 0 
Ty = -T sin y cos @ (8) 
T, = Fl cos 


Summing forces in the three component directions, one obtains: 


-T sin » sin 0 + Dp sin ¢ sin y = 0 (9) 
-T sin ~ cos 6 + Da sin ¢ cos y = 0O (10) 
-T cos yp + Dp cos 8 + (B- W) = 0O (11) 
From Equations (9) and (10): 
tan © = ‘tawy (125 
therefore: 
pee TY (13) 


(which was to be expected because the horizontal components of the drag force, 
resulting from the horizontal velocity components, determine the direction of 
the tension force at the buoy). 


From Equations 10, 11 and 13: 


a ely sin $ 
Ye seen (B-W) /D, + cos 5 Soy 


and: 


= sin 
: Ds sin (15) 


Equations 3, 4, 5, 13, 14 and 15 are solved to provide values of 0, y and T 
at the buoy. These, together with coordinate values (X, Y, Z), provide the 
initial conditions at the upper end of the mooring cable. 


B. The Cable 
Assumptions made in modelling the mooring cable are as follows: 
l. The cable is flexible and extensible. 


2. The cable diameter is not significantly reduced when stretching 
occurs, i.e. radial strain is neglected. This has implications 
for the drag cross-sectional area and for the displaced volume, 
but the effect is minute, for steel cables at least. 


3. Cable drag can be separated into normal and longitudinal 
components. 


4. No strumming occurs explicitly, although its effect can, 
perhaps, be simulated by means of a larger-than-normal drag 
coefficient. 


The solution of the cable free-body equations is more easily 
accomplished in a system of cable coordinates, rather than in the principal 
reference frame. This system is defined by first rotating the XY-plane 
(Fig. 1), through the angle 6 about the Z-axis, to bring the Y-axis into 
coincidence with the tangent to the projection of the cable on the XY plane. 
One now has a new set of axes x', y' and z', say, where: 


x" =s X cos 6 = Y sin 0 


: X sin 6 + Y cos 6 (16) 


Ke 
i] 


Next, rotate the y'z'-plane through the angle about the x'-axis to bring 
the z'-axis into coincidence with the tangent to the projection of the cable 
on the y'z'-plane. Now one has the cable coordinate axes x, y and z, where: 


x = X cos 0 = ¥Y sim ¢@ 
y = X sin 6 cos ) + Y cos 6 cos py - Z sin p Cig 
Zz = X sin ¢ sin’ ty cos Gosin ) + 2 cos -u 


This is more readily expressed in matrix form as: 


GC) AT (18) 
or: 


glLikangite (19) 


where C is the column matrix of any vector quantity in the cable coordinate 
system and I is the column matrix of the corresponding vector quantity in 
the principal reference frame. The transformation matrix A and the inverse 


matrix ee are given by: 


A = cos 6 -sin 6 0 
sin 8 cos yp cos 8 cos wp -sin wv (20) 
sin 6 sin wp cos 6 sin jp cos 
and: 
hae = cos 8 sin 8 cos yp sin, Ou sin vy 
-sin 6 cos 8 cos jp cos 8 sin p (21) 
6) -sin, cos Ww 


The matrix relationships can be used to obtain the cable system 
velocity components, namely: 


u = U cos 6 - V sin 8 
v = U sin 0 cos ») + V cos 8 cos wb - W sin w (22) 
w = U sin 0 sin ») + V cos 8 sin ) + W cos yp 


The magnitude of the resultant velocity is: 


ts a ef +: = Ro (23) 


and its orientation in the cable system is defined by the angles n in the 
xz-plane and in the xy-plane. The following relationships hold: 


- u Vv W ise vad 
R _ SSS __—=— = iS ee 
sin n sin A sin n cos xX cos n sin n (24) 
so: 
4 TES ig 
Nn = cos f S. C25) 
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(26) 


An element of cable length dS is considered. This element lies in 
the direction of the z-axis and its location in the principal reference frame 
can be found by integrating along the cable length S over the directional 
derivatives: 


OX wine 4 ‘ : 

a8 sin 6 sin wv 

aye *.” , 

qs cos 6 sin w C27) 
WZ. as 

GC cos 


In considering the hydrodynamic drag forces on the cable, the 
Independence Principle of fluid dynamics is assumed to apply (Schlichting, 
1960), i.e. the normal drag force depends only on the normal velocity (the 
component in the xy-plane). This force comprises both pressure drag and skin 
friction. The longitudinal drag force consists solely of skin friction, with 
the stress proportional to the longitudinal velocity. Thus, two independent 
forces exist: 

1: 
2 
f = (f. +f a o (a + hy 
(28) 


If q, r and s are drag forces per unit length in the x, y and z directions, 
respectively, then: 


f = q dS 
Pe |. Eads 29 
9 (29) 
£ = gs dS 


To solve for q, r and s, the usual convention is followed of assigning a non- 
dimensional drag coefficient as the proportionality constant in Equations 28. 
Thus: 


2 poh 7A 
fy = Ch 5 (tie aaa a) Leds (30) 


where D is the cable diameter and Ch incorporates the effects of both pressure 


drag and normal skin friction. Cy can be written as: 
C.F fuze ves (31) 


where the pressure coefficient C, is measured at n = 90°, and the skin 


P 
friction coefficient CoN is a function of Reynolds number R.: Also: 
fe Pee Roe dS (32) 
is FL 2 


where the longitudinal skin friction coefficient is a function of both the 
angle n and Ro» as: 


1 1 
G = “k Saag cos n sin” n (33) 


The constant k is determined empirically. The appropriate form for skin 
friction coefficients is open to some conjecture. For a fuller discussion 

of cable drag, the reader is referred to Bell (1977). Now the drag components 
per unit length can be determined from the foregoing equations, as follows: 


in A 
q ds ds D2 y 
£ fo COSA 
CT iy BELO NOU ERS gh [ek Sa ae a 
r 4s ds Cy 5 DR sine ae cos x (34) 
f be 
ees te Sele Hae? i 2 
de ds Sr Rivie alley 


Gravity forces on the cable consist of a weight force w per unit 
length and a buoyancy force b per unit length. The combination acts in the 
Z-direction, with a magnitude of (b-w)dS on each elemental length of cable. 
Resolving this in the cable coordinate system: 


Seen 
By = (w-b) sin wy dS (35) 
g = -(w-b) cos wp dS 


Tension components at the upper end of the cable are, in the 
principal reference frame: 


Ty = (T+dT) sin(\yt+dy) sin(6+d6) 
Ty = (T+dT) sin(ytdy) cos(6+d6) (36) 
Ty = (T+dT) cos( +d) 


Transforming to the cable coordinate system by use of the previous matrix 
relationships, and ignoring second order quantities, one obtains: 


toe = T dé sin wp 
38s = T dy (37) 
= T+ 
tue T+dT 
At the lower end of the cable: 
tix = Dye = 0 
(38) 
t = -T 


Summing all of the foregoing components in each of the cable 
coordinate directions gives a set of three equations: 
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Wis 2) 4 Verges 
ds T sin w 


“ = ~3 {r + (w-b) sin y} (39) 
= = -s + (w-b) cos yp 


These three equations, together with the defining equations for the drag 
force components, and with Equations 27, are solved simultaneously to obtain 
all of the required cable parameters. 


An inextensible cable has been assumed in the foregoing. If an 
elastic cable is desired, then it is necessary to make the following 
substitutions, where e€ is the elongation per unit length: 


dS + dS (l1+e) 
(40) 
(w-b) > (w-b)/(1+e) 


Here, the implication is that the weight of one foot of stretched cable is 
w/(l1+e) so there is no change in total weight. Further, the buoyancy per 
foot of stretched cable happens also to be taken as b/(l1+e), so that some 
compensation for the effect of radial strain on the displaced water volume is 
included, although not in the correct formal manner. Any error in this term 
is small compared to the other forces involved, so it was not considered 
necessary to formalize it. Also, the effect of radial strain on the drag 
force term is neglected, again because of its relative insignificance. To 
approximate the elastic characteristics of various cable materials, it is 
convenient to use a second-order curve: 


kb +6005 a. a + ave (41) 


where T is the tension and Ags ay» ay are constants to be determined from 


the material specifications supplied by the cable manufacturers. An expanded 
discussion of this subject is available in Bell (1977). The set of equations 
for a stretchable cable is, then: 
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clea <q/(i¥Fe) 


ds T sin yp 
& = -z {r(1lt+te) + (w-b) sin yp} 
< = -s(lte) + (w-b) cos yp 
(42) 
& = alae os ingepis inal 
= = =(1+6)" cos“8- sin 
& =) Clbhe) cose! 


C. The Instruments 


Instrument packages of various shapes and sizes may be attached to 
the cable at various locations along its length. Here, as for the buoy, only 
the simplest case is considered. All of the forces act through the centre 
of gravity and the drag coefficient is not a function of aspect. The 
instrument has a weight Wy in air, a buoyancy Bre a drag coefficient Cor 
and a cross-sectional area Als Its presence in the system causes incremental 
changes, at the point of attachment, in line angle, line azimuth and tension, 
given by Ay, A® and AT, respectively. The drag on the instrument is given 
by: 


Déct=Ene R|R| (43) 


2 
I DT o°"t 


In the principal reference frame, the drag force components are: 


Diy = Dy sin ¢ sin y 
Dy = D, sin ¢ sin y (44) 
D aE TL COS. 


IZ ii 
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The tension components at the upper cable attachment point are: 


Toy = T sin wv) sin 6 
Toy = T sin y cos 6 (45) 
Toyz = T cos 0 


At the lower cable attachment point, they are: 


Thy = -(T-AT) sin(vt+Ay) sin(6+A6) 
Ty = -(T-AT) sin(ytAy) cos(6+A8) (46) 
Tz = -(T-AT) cos(vtAp) 


The gravity force is: 


(47) 
He. do any ee 
Summing the components in the three coordinate directions, and 
solving the resulting equations, gives the new cable azimuth, angle and 
tension due to the presence of the instrument as: 
D + °T sin sin 0 
a -l IX 
pa cag al D + T sin wv cos 6 Os) 
LY 
rear ean) Diy + T sin Wy cos 8 Le 
cos(6+A8) {D_. + T cos »y + BL - W_} 
We 1 I 
D.~-+ T cos ) *B_7—w 
T-AT = ee : i (50) 


cos (yt+Ay) 


Some instrument packages may be very buoyant or very heavy, so the simplifying 
approximations for small changes in angle cannot be made here. 
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D. Method of Solution of System Equations 


The equations for the complete buoy-cable system are first-order 
non-linear and form a boundary-value problem requiring boundary conditions at 
each end of the cable. The buoy boundary condition is satisfied by assuming 
an elevation for the buoy. This, in turn, fixes the velocity and permits 
calculation of the conditions at the upper end of the cable. The anchor 
boundary condition is that the elevation of the anchor must coincide with that 
of the sea bottom, within acceptable limits. The governing equations are 
solved for the unknown coordinates, angles and tension, by a step-wise 
integration along the cable, beginning at the buoy end, to a point of 
discontinuity resulting from the concentrated loads caused by the presence of 
an instrument package. The effect of the forces on the package is added to 
the cable forces and the integration is continued in a similar fashion until 
the lower end of the cable is reached. If the anchor boundary condition is 
not satisfied within the specified error limits, then a new buoy position is 
chosen and the procedure is repeated. 


Computer Solution of the Buoy-Cable System 
A. The Buoy 


No specific geometry having been used in formulating the model, the 
only buoy-related input data required by the computer program is the weight 
in air of the buoy, the gross buoyancy (i.e. the weight of water displaced) 
and the product of the drag coefficient and the cross-sectional area presented 
to the flow. If a characteristic dimension is also included in the input 
data, then the Reynolds number will be calculated. An initial elevation is 
assigned to the buoy to fix its position with respect to the velocity field 
and to the bottom. This elevation is assumed to coincide with the bottom of 
the buoy, which is the point of attachment of the cable. Calculation of the 
boundary conditions for this end of the cable can then be carried out 
directly. Here, and elsewhere throughout the program, numerous tests are 
incorporated to avoid divide-by-zero situations and similar problems. 


B. The Instrument 


‘The program permits various kinds of instrument packages, including 
floats for additional buoyancy, to be inserted in the system. However, they 
are represented simply as point forces and should be positioned at multiples 
of the basic cable segment length used for integration along the cable 
(discussed in the next section). The instrument locations, in length units 
along the cable from the buoy, are read into the computer, converted into an 
integer number of cable segments and stored in an array capable of holding 
fifteen different values. Because of the conversion, the instrument locations 
are still correct when the cable is stretched. Other data which must be 
provided are the weight in air, the gross buoyancy, the cross-sectional area 
presented to the flow and the drag coefficient. Any shackles or other fit- 
tings located near an instrument should have their weights accounted for by 
including them in the instrument weight. Before each integration is made, 
the program checks to see if there is an instrument at the upper end of the 
cable segment. If so, the velocity is determined and the instrument drag 
calculated. Then the resulting increments in cable tension and angles are 
combined with the preceding values of these variables. 
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C. The Cable 


The integration along the cable is carried out over a series of 
cable segments which are of equal length when unstretched, this length being 
a sub-multiple of the total cable length. The choice of segment size depends 
on the problem parameters, such as the complexity of the velocity profile, 
the number and location of instrument packages, etc. The integration routine 
uses the Runge-Kutta method with error control. An external subroutine 
specifies the functions required for the evaluation of the derivatives given 
in Equations 42. The integration routine also requires the specification of 
an error tolerance on the integral, a step size, and a value of the 
independent variable at the end point of the integration. Double precision 
is used throughout to reduce the possibility of cumulative rounding errors. 


For a stretchable cable and step-wise integration, the end point 
must correspond to the strained length of each cable segment. This requires 
that the independent variable be stretched as well, i.e. Equation 42 is 
divided through by the quantity (l+te), giving (l+e)dS as the independent 
variable instead of dS. In the program, the strained unit length (lt+e) is 
represented by the function STRCH. It is calculated on the basis of the 
tension at the upper end of a cable segment. For greater accuracy, this 
tension could include one-half of the tension increment for the previous 
segment. Generally speaking, however, the increment is usually no more than 
a few tenths of one percent of the tension, so this correction is ignored 
here. If required, the accuracy can be improved by reducing the segment 
length. Any necessity for this is easily checked by running the model for 
two different segment lengths. 


The boundary conditions for the first cable segment are provided by 
the buoy calculations, and for subsequent segments by the results of the 
integration over the preceding cable segment. The velocity used for calculat- 
ing cable drag is that which occurs at the midpoint of the segment, the mid- 
point elevation being estimated on the basis of the cable angle for the 
preceding segment. An iteration on the cable angle would result in a slightly 
more accurate midpoint location and, therefore, a slightly more accurate 
velocity (only when the velocity profile is not uniform). This was not deemed 
worthwhile since, for the calculation of drag, the midpoint velocity is 
assumed to apply over the whole segment in any case. One could incorporate a 
varying velocity (when such is the case) throughout the integration over a 
cable segment, but at the expense of a considerable increase in program 
complexity. Again, if concern is felt about the outcome of the procedure as 
presently used, especially when the velocity is changing rapidly with depth, 
the program can be run with small segment lengths. Some trial runs in which 
the segment lengths were varied by factors of two or more, all other 
variables remaining unchanged, showed no variation in the results. Thus, the 
indication is that the simple procedure used to obtain the velocity for the 
drag determination is adequate. Also, while only an average velocity value 
is used in the drag calculation, the dependence of the drag on the actual 
cable angle is included throughout the integration along a cable segment since 
the angle is one of the dependent variables. 


The mooring line used in the model is not restricted to one size or 
one kind of material throughout its length. Presently, ten different cable 
types can be incorporated. Information on these, including the distance from 
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the buoy of the lower end of each type, is read into an array. Then (as is 
done for the instruments), the end position is converted into an integer 
number of cable segments, the conversion permitting an easy check on type 
changes as well as allowing cable stretch without causing any additional 
problems. The form of the longitudinal drag equation used with each type is 
optional, the choice being either that for smooth cable, for wire rope, or 
no axial drag at all. The selection is made by assigning an appropriate 
logical flag in the input data. 


D. The Velocity Profile 


Provision is made in the model for incorporating a three-dimensional 
velocity field around the cable in terms of the components U, V and W. The 
speed versus depth relationships for the components are represented by first- 
order profile segments of the form: 


U = ay 5 by Z 
Ve= a, + b, Z (51) 
W = a, a b, ZL 


where Z is the elevation above the sea floor. Up to ten such segments can be 
pieced together for each component, permitting the simulation of a fairly 
complex velocity field. The profile segments can be any convenient size and 
do not have to correspond to multiples of the cable segment length. In the 
program, the velocity profile information is stored in an array, the first 
element in each row being the lower depth limit of applicability (measured 
from the sea floor) of the three-dimensional segment whose coefficients 

(a)> by> aos bos ax, b,) are given by the other elements of that row. A 


conditional branch on the depth is used in the various calculations to obtain 
the correct velocity at a given depth. 


E. The Iteration Procedure 


The solution of the buoy-cable model is based on making an estimate 
of the buoy position with respect to the bottom and, proceeding from there, 
integrating along the cable to the anchor. If the anchor elevation doesn't 
coincide with the bottom elevation, within specified limits, the estimate of 
the buoy position is modified and another integration is performed. A simple 
additive modification of the buoy position is made for the first iteration, 
using the error in the anchor location. Subsequent iterations use a 
procedure which is, perhaps, more complex than necessary for subsurface 
moorings since it is a carry-over from the two-dimensional model, where the 
complexity was required for truculent cases involving surface floats. It 
does seem to provide a more rapid convergence to the correct solution than 
does any simpler scheme tried. The program halts when more than fifteen 
iterations are required, but subsurface moorings seldom take more than one 
or two. 
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F. Program Output 


The first portion of the program output lists information 
associated with the integration and iteration routines. In particular, the 
various estimates of buoy elevation and the resulting anchor position are 
given for each iteration. Following this there is a listing of all pertinent 
input and output data about the buoy-cable system and the velocity profile, 
including items such as instrument package drag and position coordinates. 


Next, a printer plot of the coordinate system is given as in atde- 
memotre when perusing the various angular quantities in the tabulation which 
follows. This table gives the coordinates, tensions, forces and velocities 
corresponding to the solution of the system equations for each cable segment. 
The cable drag values given in this table are for unit cable length. The 
drag value listed in the j-th row applies to the segment whose lower-end 
coordinates are also listed in the j-th row. 


G.v Verification 


The model has been exercised in all eight octants with the same 
contrived three-dimensional data set and appropriate symmetrical results 
obtained. This implies that there were no basic programming errors in 
implementing the three-dimensional equations since, if errors had been made, 
one might expect some asymmetries to appear. The model has also been applied 
to field data previously used to verify the two-dimensional program and 
identical results obtained. Taken together, the results of these tests have 
inspired confidence in the use of the model for three-dimensional predictions 
although, as yet, no adequate field data set has been available for a proper 
three-dimensional verification. In particular, adequate velocity profile 
information is required because of the square law dependence of drag on flow 
speed. 


Conclusion 


A model of a static single-point subsurface mooring system has been 
developed for predicting system geometry in a three-dimensional velocity 
field. It incorporates a flexible, extensible cable whose characteristics 
may be varied over different portions of the system and it permits the 
inclusion of instrument packages at various locations along the cable. The 
accuracy of results obtained from the model depends in large measure on the 
provision of accurate velocity information. 
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Appendix A_- Sample Computation 


The program is here applied to a system comprising a subsurface 
float, two instrument packages and a 100 ft length of cable which is mainly 
wire rope with a short segment of nylon line at the lower end. The FORTRAN 
coding form for this example is given in Table I. 


The first data card provides the buoy specifications, with ten 
columns allotted for each of the following: the gross buoyancy (weight of 
water displaced - 600 lbs in this instance), the buoy weight in air (300 lbs), 
the diameter (3 ft) or some other characteristic length, and the product of 
drag coefficient and cross-sectional area presented to the flow (4.2 sq ft). 


The next cards in sequence provide information on the cable. Here, 
two different types are used. The second data card is for '' diameter wire 
rope whose lower end is at a distance of 90 ft along the cable from the buoy, 
as shown in Column 1-7. The next five sets of seven columns each are used to 
specify, in the order given, cable gross buoyancy per unit length 
(0.022 lbs/ft), weight in air per unit length (0.100 lbs/ft), diameter 
(0.021 ft), drag coefficient (1.4) for flow in a direction normal to the 
cable and drag coefficient (0.0) for flow parallel to the cable axis. Begin- 
ning in Column 43, there are three sets of ten columns each which can contain 
the coefficients for a second-order fit to the fractional stress-strain 
relationship for the cable. In this instance, since wire rope is being used, 
Hooke's Law is obeyed so only first-order coefficients are required (1.0, 


oea xX ney 0.0). Columns 73-75 contain logical variables used for selecting 


the form of the longitudinal drag equation. If a "T" appears in Column 73, 

a form appropriate to smooth wire is used; if in Column 74, the form is for 
rough wire; if in Column 75, longitudinal drag is ignored altogether (as 
here). The third data card represents a segment of }5"" diameter nylon line 
whose lower end is 100 ft from the buoy. Buoyancy, weight, diameter and drag 
coefficients are, respectively: 0.045 lbs/ft, 0.050 lbs/ft, 0.036 ft, 1.4, 

8 


> 


0.1. The stress-strain coefficients are 1.045, 1.125 x 1074 and 1.389 x 10. 
as determined from a graph supplied by the manufacturer of the line. The "T"' 
in Column 74 shows that the rough wire longitudinal drag equation is being 
used for this case. Up to ten different cable types may be represented in the 
model, each with its own data input card. If the number of types used is less 
than ten, the sequence must be terminated with an end-of-file indication. 
Hence, the fourth card here contains @ EOF. 


Next, the instrument package data is supplied. The fifth and sixth 
cards contain, in Columns 1-10, the locations of the instruments in distance 
along the cable from the buoy. Instruments must be positioned at multiples 
of the basic segment length, as mentioned in the text. In this example, the 
length is 5 ft and the instruments are at 85 and 90 ft. Then Columns 11-20 
of each card contain the gross buoyancy (50 lbs and 20 lbs, respectively, for 
the two instruments), Columns 21-30 show the weight in air (15 and 50 tbs), 
Columns 31-40 hold the cross-sectional area for the drag calculation (0.8 and 
1.5 sq ft) and Columns 41-50 contain the drag coefficient (0.6 and 6 
Since the number of instrument data cards is less than the maximum number of 
fifteen, the next card must be an end-of-file indicator. 
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Following on, another five cards are used to provide the velocity 
profile information. Each card specifies a particular segment of the profile 
by giving, firstly, the lower limit of applicability of the segment in 
distance above bottom and, secondly, the coefficients for a linear velocity 
vs depth equation for each component. Thus, the eighth card indicates a 
uniform velocity, at elevations in the water column greater than 60 ft above 
bottom, having a U-component of 2 ft/sec, a V-component of 5 ft/sec and no 
W-component. The ninth card shows that, between 50 and 60 ft above bottom, 
the velocity components are U = -13 + 0.25 Z, V = 5, W = 0, and similarly for 
cards ten through twelve. The easiest way of determining the coefficients is 
to graph each component of the desired velocity profile, for clarity, then 
obtain a zero or first-order fit to each segment separately. Ten velocity 
profile segments are provided for in the program; if fewer than this are 
used, an end-of-file indication must be inserted. 


The fourteenth, and final, data card contains the desired number of 
cable segments (20) in Columns 1-10, the cable length (100 ft) in Columns 11- 
20 and the initial guess for the buoy elevation (100 ft) in Columns 21-30. 
The initial guess can be any reasonable value because the iteration procedure 
usually ensures a rapid closure on the correct position. The permissible 
iteration error, which is the position accuracy for the anchor end of the 
cable, is provided in Columns 31-40. Here it is 0.25 ft. The integration 
subroutine requires the specification of an integration error tolerance, 
given here as 0.01 in Columns 41-50. Generally, a value of 0.01 to 0.1 has 
been satisfactory in all of the modelling done to date. Finally, Column 51 
contains a logical variable which, if given as "'T", will result in the 
printing of various data as a diagnostic aid, when required. 


When the program is run, various listings are produced, as shown 
in Table II for the present example. These listings are self-explanatory. 


TABLE II. 


NO. 


Husa 


OUTPUT LISTING FOR MOORING PROGRAM EXAMPLE 


OF CABLE SEGMENTS = 


INTEGRATION ERROR TOLERANCE 


PERMISSIBLE ITERATION ERROR 


INITIAL GUESS FOR BUOY POSITION = 


ZED 


100.00 
89-76 


89.76 
90.38 


90.38 


NO. 


VELOCITY PROFILE 


LOWER LIMIT 
(FI ABOVE BOTTOM) 

o0e 

S06 

40-6 

256 

O- 


Y(6) 


10.24 
10.24 


-.66 
-.66 


00 


2.00 
-13.00 
-13.00 

-3.00 
-3.00 


PIPES SS 


OF ITERATIONS REQUIRED WAS 2. 


20. 
e01000 

025 FT 

100.00 
ZMAX ' ZLOW 
100.0u 79.51 
100.0U 89.76 
e000U0¥Z 5.00 + 
°©25000«Z 5.00 + 
250002 -9.00 + 
e000U0 «2 -9.00 + 
e000U0*«2Z -2.00 + 


FT ABOVE ROTTOM 


YMAX 


10.24 


10.24 


000002 
00000*2 
-28000*2 
e28000*2Z 
e00000*Z 


00 
00 
00 
00 
200 


= 


+++ 4+ 


YMIN 


-10.24 


- 266 


000002 
e00000*«2Z 
e00000*2 
-00000*«2 
e00000«2 
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TABLE II. (continued) 
BULY PARAMETERS : 
GRUSS BUOYANCY = 600.000 LBS 
WEIGHT IN AIR = 300.000 LBS 
DIAMETER = 3.000, Fl 
CRUSS=SECTIONAL AREA X DRAG COEFFICIENT = 4.2005 150 ET 
BUOY COURDINATES ( 14.50 » G2N1T) 90.38 ). 
VERTICAL EXCURSION = 9.62 FT 
BUY DRAG = 121-80 LBS 


CABLE PARAMETEKS = TYPE 1 3 


GROSS BUOYANCY = °0220 LBS PER FT 
WELGHT IN AIR = «1000 LBS PER FT 
DIAMETER = ©0210 Fas 

NORMAL DRAG COEFFICIENT = 1.400u 
AXLAL DRAG COEFFICIENT = ©0000 


ELASTIC COEFFICIENTS: AO= .1000+001 + Al= .2400-005 ; A2= .0000 
LENGTH = 90.000 FT 


CABLE PARAMETERS = TYPE 2 ; 


GRUSS BUOYANCY = e0450 Efsisy detest Ley 
WEIGHT IN AIR = ©0500 EBSSRERSE Ti 
DIAMETER = 20360 FT 

NORMAL DRAG COEFFICIENT = 1.400u 
AXLAL URAG COEFFICIENT = -1000 


ELASTIC COEFFICIENTS: AV= 610454001 + AL= 1125-003 + A2=-.1389-007 
LENGTH = 1U.0U00 Fall 


TOTAL LENGTH OF CABLE = 100.000 FTe 


INSTRUMENT PARAMETERS = NO. 1 3 


VISTANCE ALONG UNSTRETCHED CABLE FROM BUOY = BSe FT 
BULVYANCY = 50.000 LBS 

WEIGHT IN AIR = 15.000 LBS 

CRUSS=SECTIONAL AREA = -800 SQ FT 

DRAG COEFFICIENT = e600 

INSTRUMENT ORAG IS 6624 LBS 

INSTRUMENT COORDINATES ( 1.08 » 6633 » 14.44 ),. 
VERTICAL EXCURSION = Alois) [PY 


INSTRUMENT PARAMETERS = NOew 2 3 


OISTANCe ALONG UNSTRETCHED CABLE FRKOM BUOY = 90. FT 
BUOYANCY = <0.000 LBS 

wELGHT IN AIR = 50.000 LBS 

CROSS=SECTIONAL AREA = 1.500 SOeET 

DRAG COEFFICIENT = 1.1uU0 

INSTRUMENT DRAG IS 21.45 LBS 

INSTRUMENT COORDINATES ( ©5355 » 4.32 »¢ 9.89 ). 


VERTICAL EXCURSION = eek EF 
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TABLE II. (continued) 


AN ATTEMPT TO ILLUSTRATE THE COORDINATE SYSTEM USED FOR THE 3=D MOORING MODEL 

IS PRESENTED GELOWe THE ANGLES ASSOCIATED WITH THE CABLE VECTOR ARE *ANGLE* & *AZIMUTH*> 
AND WITH THE CURRENT VECTOR ARE *TILT*® & *DIRECTION*®. NORTH LIES IN THE DIRECTION 

OF THE Y-AXIS AND EAST IN THE DIRECTION OF THE X=AXIS. 


Zz 
& 
W 
* 
* 
x Le 
* B 
* A ° 
* ANGLE Cc 
*efo 0 e e 
ot . E 
* AoL : 
“p* B 
TILT . * A ° 
(e ° * C 
a ° x é 
R > x (S 
r R * ie ry 
E * B 
é N * A ° 
Liat Ks 
** & ee & * * ke eK kK KOK Ke & Y 
ry 2 4 * 4 é V 
r * 6 a e e r 
* e e 
e * ° e @ AZIMUTH e 
* We e ° ° 
e * ° DIRECTION 7.0 e 
* 
e * e e e 
* 
* 6 . e e 
U 
& e e e e 


(continued) 
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Appendix B_ - Program Source Listing 


A complete listing for the main program and subroutines is given in 
the following pages. 
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